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Preface 



This is a book on empirical number theory concentrating on the analysis of 
number sequences. Its focus is on a small part of a very large number of 
integer sequences defined by Florentin Smarandache. The author has, 
however, when appropriate included some other of his research results which 
organically belongs to this area. The content is organized into chapters 
according to the main considerations to make when programming the 
analysis. They are not mutually exclusive. 

As in my previous book Surfing on the ocean of numbers - a few 
Smarandache Notions and Similar Topics an attempt has been made to 
present results so that they are easy to understand and not to burdensome to 
read. There are many tables, some of which may be used for reference, but in 
general they are there to show the overall results obtained at a glance. In most 
cases it is the way in which sequences behave, not the individual figures, 
which is of importance. Some graphs have been included to illustrate 
important findings. 

Some of the results in this book were presented by the author at the First 
International Conference on Smarandache Notions in Number Theory , 
August 21-23, 1997, Craiova, Romania. 

References have been given after some of the chapters. However, constant use 
has been made of the following Smarandache source materials: 

Only Problems, Not Solutions!, Florentin Smarandache. 

Some Notions and Questions in Number Theory, C. Dumitrescu, V. Seleacu. 
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Illustrations, graphics, layout and final editing up to came ra ready form has 
been done by the author. Most tables have been created by direct transfer from 
computer files established at the time of computing to the manuscript so as to 
avoid typing errors. The results often involve very large numbers which are 
difficult to accommodate in pie-designed formats. It has therefore often been 
necessary to use several lines to represent a number in tables as well as text 
All calulations have been carried out using UBASIC, ver. 8.87. 

I am grateful to Dr. R. Muller who has given all possible help and 
encouragement during the work on this book. He and his colleagues at the 
American Research Press have at all times facilitated the work through rapid 
e-mail communications. 



Finally, it is thanks to the patience and understanding of my wife Anne-Marie 
that this book has come about But the winter months have passed, the book is 
finishe d and I am no longer going to be lost among manuscript pages. 
Summer is around the comer and the Swedish nature is waiting for us. 



Paris, April 1998 
Henry Ibstedt 
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Chapter I 

Partition Sequences 



LI Introduction 

We will study the partition of the sequence of digits 

of positive integers into two or more groups of digits whose corresponding 
integers may or may not possess a certain property. 

Example: 

17 | 7 | 101 

is a partition of the 6-digit integer 177101 into three groups of digits (parts) 
each of which represents a prime, whereas 
1 | 77 | 10 | 1 

is a partition of the same integer into four parts none of which is a prime. 

Obviously there are many ways in which and integer can be partitioned into 
groups of digits. The kind of partition described above may also be referred to 
as a de-concatenation in order to distinguish it from the classical concept of 
partition of an integer as first developed by Euler and later illustrated by 
Ferrer [1]. However, the kind of partition we are dealing with will be clear 
from the context. We will now use the classical partition of an integer as a 
means to help us formulate a strategy for partitioning sequences of integers in 
the sense defined above. 

The partitions of 4 are: 1+1+1+ 1=2+1 +l=2+2=3+ 1=4 with the corresponding 
Ferrer diagram: 



□□□□ □□□ □□ □□ □ 

□ □□ □ □ 

□ □ 

□ 

1 2 3 4 5 

Diagram L The Ferrer diagram for 4. 
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The last partition corresponds to 4 itself which is of no interest to us. 
Furthermore, we will link the vertical representations to form the pattern 
shown in diagram 2: 

□□□□ □□□ □□ czzo 

12 3 4 

Diagram 2. The Ferrer diagram for 4 rearranged and with partition numero 5 
excluded. 

In our applications we will also need to consider the order in which the 
partition elements occur. Diagram 2 shows the seven different ordered 
partitions possible for 4 and how they form a partition pattern for the 
partition of the 4-digit integer 9164 into groups of integers. 

0000 cues 000 0bc m [mu] mnm mn*n 

1 2 3 4 5 6 7 



Diagram 3. The partition pattern fora 4-digit integer. 

From this we see that the integer 9164 can be partitioned into perfect squares 
in exactly two ways, 9, 16, 4 (partition numero 3) and 9, 1, 64 (partition 
numero 4). This is how the mapping of and integer onto a partition pattern 
helps us study the properties of each partition element 

Only partitions 2 and 4 in diagram 2 con tain partition elements of unequal 
sizes. These give rise to an increase of the number partition patterns through 
arrangement of the partition elements. It is easily understood that in a 
partition erf 8 we have 2 parts of the type L— J and 4 parts of the type □ 

then the number of arrangements will be * 5 * The integer 8 has 21 

partitions of the type shown in diagram 2 giving rise to 127 ordered 
partitions. This could be calculated by considering each arrangement as we 
did above. However, we will soon arrive at this in a different way. 



s 




The above considerations will be useful in order to express the partitions of 
integers in explicit form. When we consider the partition of a very large 
number of integers to find out how many of them possess a certain property 
we will be able to use a different approach. 

Let’s represent an n-digit integer t in the form t=ai+a 2 * 10+a 3 *10 2 + ... a n -10 n ‘ 1 . 
We introduce the following definitions and computer analysis related 
concepts in UBASIC. 

5(t)=l if t possesses a certain property otherwise 5(t)=0. 

\ symbolizes integer division, example 23678\100=23 

res symbolizes the remainder of the last performed integer 
division, in the example above res=678. 

We will now take a look at the question of leading zeros. Is 3 | 07 a partition 
of 307 into two primes? Is 8 | 027 a partition of 8027 into two cubes? The 
author prefers to have a unique representation of an integer, 7 is a 1 -digit 
integer written 7 not 07 and 3 3 =27 is written 27 - not 027. The computer, 
however, interprets 07 as 7 and 027 as 27. To avoid integers with leading 
zeros the function e(j,r) is defined as follows. 

if KlQ^ 1 otherwise s(j,r)=l, where r and j are context 
depending variables. 

We will now use these functions to formulate an algorithm for calculation of 
the number of n-digit integers which can be partitioned so that the partition 
elements have a given property (being primes, squares, cubes, etc) 

Lets denote the number of partitions of t by p(t). 

1- digit integers : l^t<9. Obviously there is nothing to partition, p(t)=0. 

2- digit integers : 10<t^99. p(t) can only assume the values 0 and 1 and is 
given by 

p(t)=5(t\10)-5(res) 
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3-digit integers : 100£t£999. In the algorithm below use is marie of the fact 
that we know p(t) for t<100. 

p(t)=S(t\100)-E< 1 ,res)-{5(res>+p(t\100)>+5(t\10H5(res>+p(t\10)), 

where we know that p(t\10)=0. It has been included in order to present a 
general algorithm for n-digit integers. 



n-dieit integers : Hf’stSlO 0 -!. 



It- 1 

p(t) = ]£<?(/ \ Hf~ k )e(n -k- \,res)(S(res) + p(res)) 



k = 1 



( 1 ) 



We can now use (1) to calculate the number of partition patterns for n-digit 
integers. To do this we pot 8(t\tf**)=l (the property of being a partition 
element), e(n-k-l,ies)=l (a partition element may begin with a zero) and 
p(res)=p(k). This results in 



B-l 

Pin) = X( 1 + Pi *)) . whcre P(1)=0 

k=\ 

Evaluation of the recursion formula gives p(n)=2 D " 1 -l resulting in the 
following taMe of partition patterns. 

Table 1 . Partition patterns. N = number of patterns. 



n-digit integer 1 2 3 4 5 6 7 8 

N 0 1 3 7 15 31 63 127 



The algorithm has been implemented in UBASIC on a Pentium 100 Mhz 
portable computer to examine various partition sequences. 



to 





L2 The Smarandache Prime-Partial Digital Sequence 



Definition: The prime-partial-digital sequence is the sequence of prime 
numbers with the property that each member of the sequence can be 
partitioned into groups of digits such that each group is a prime. 

Table 2 shows the first one hundred terms of this sequence. 

We will see in Chapter H4 that the Smarandache Prime-Digital Sub- 
Sequence is infinite. The Prime-Partial Digital Sequence is in fact a subset of 
the Prime-Partial Digital Sub-Sequence and is therefore a fortiori infinite 
which settles affirmatively a conjecture by Sylvester Smith. [2]. We will make 
observations on the number of ways in which prime integers can be 
partitioned into primes. How close is the maximum number of ways this can 
be done to the theoretical maximum given in table 1 and what do these 
partitions look like? 



Table 2. The first one hundred members of the prime-partiakigit sequence. 



23 


37 


53 


73 


113 


137 


173 


193 


197 


211 


223 


227 


229 


233 


241 


257 


271 


277 


283 


293 


311 


313 


317 


331 


337 


347 


353 


359 


367 


373 


379 


383 


389 


397 


433 


523 


541 


547 


557 


571 


577 


593 


613 


617 


673 


677 


719 


727 


733 


743 


757 


761 


773 


797 


977 


1013 


1033 


1093 


1097 


1117 


1123 


1129 


1153 


1171 


1277 


1319 


1327 


1361 


1367 


1373 


1493 


1637 


1723 


1733 


1741 


1747 


1753 


1759 


1777 


1783 


1789 


1913 


1931 


1933 


1973 


1979 


1993 


1997 


2113 


2131 


2137 


2179 


2213 


2237 


2239 


2243 


2251 


2267 


2269 


2273 



The property that we will examine with the function 5(t) is whether t is a 
prime or not. In UB ASIC this is done in the following program. 

5(t): 10 ,D(t) \D(t) is equivalent to 8(t) 

20 if t=0 then z=0 :goto 40 

30 if nxtprm(t-l)=t then z=l else z=0 

40 retum(z) ‘z is the value of D(t) 
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The function e(j,r) takes the following form in UBASIC. 

e(j,r): 10 .L(j%,r) *JL(j%,r) is equivalent to 

e(j,r) 

20 if KlO A (j%-l) then y=0 else y=l 

30 retura(y) ‘y is the value of e(j,r) 



These two functions woe made work together with the algorithm (1) in a 
UBASIC program to produce the results shown in tables 3 and 4. 

Table 3. ratio m/q where m is the number of rvcSgit primes which can be 
^ ffcsIBfoned into primes and q is the number of n-digif primes. 



n=2 


n=3 


n=4 


n=5 


n=6 


n*7 


n~8 


“U 4 ? 


51 


383 


3319 


27111 


234229 


2040170 


q 2^ ; 


143 


1061 


8363 


68906 


586081 


5096876 


m/q 0. 190o 0 


.3566 0.3610 ( 


} . 3969 


0 * 3934 


0.3997 


0.4003 




2 3 4 5 6 7 8 



Diagram 4. The ratio m/q for n-digit integers (n=2,3....8). 
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Table 4. The number of n-digit primes that can be partitioned into primes in p 
different ways. (The tail for n=8 is interposed.) 



p 


n=2 


n=3 


n=4 


n=5 


n=6 


n=7 


n=8 


0 


17 


92 


678 


5044 


41795 


351852 


3056706 


1 


4 


40 


251 


1986 


14789 


118379 


972095 


2 




10 


83 


769 


6306 


54102 


436009 


3 




1 


40 


309 


2811 


25892 


228428 


4 






5 


135 


1445 


14205 


139755 


5 






4 


67 


758 


7676 


82265 


6 








33 


424 


5012 


55793 


7 








15 


250 


2922 


35326 


8 








3 


131 


1877 


24401 


9 








1 


78 


1258 


17131 


10 








0 


46 


873 


12424 


11 








0 


22 


550 


8724 


12 








1 


27 


481 


6671 


13 




P 


n=8 




12 


277 


4841 


14 




38 


16 




6 


212 


3573 


15 




39 


17 




3 


147 


2781 


16 




40 


8 




0 


115 


2052 


17 




41 


13 




0 


62 


1574 


18 




42 


5 




3 


59 


1244 


19 




43 


4 






33 


992 


20 




44 


8 






34 


811 


21 




45 


14 






18 


617 


22 




46 


3 






12 


500 


23 




47 


9 






6 


417 


24 




48 


6 






10 


311 


25 




49 


2 






4 


264 


26 




50 


2 






4 


214 


27 




51 


4 






2 


154 


28 




52 


1 






3 


150 


29 




53 


0 






0 


94 


30 




54 


4 






1 


83 


31 




55 


0 






2 


85 


32 




56 


3 






0 


65 


33 




57 


0 






0 


51 


34 




58 


0 






0 


42 


35 




59 


0 






0 


51 


36 


! 


60 


1 






0 


30 


37 












1 


32 
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Table 5. All possible partitions of 23733737 into primes. 



1 


2 


3 


7 


3 


3 


7 


3 7 


2 


23 


7 


3 


3 


7 


3 


7 


3 


2 


37 


3 


3 


7 


3 


7 


4 


2 


3 


73 


3 


7 


3 


7 


5 


2 


3 


7 


3 


37 


3 


7 


6 


2 


3 


7 


3 


3 


73 


7 


7 


2 


3 


7 


3 


3 


7 


37 


8 


23 


73 


3 


7 


3 


7 




9 


23 


7 


3 


37 


3 


7 




10 


23 


7 


3 


3 


73 


7 




11 


23 


7 


3 


3 


7 


37 




12 


2 


37 


3 


37 


3 


7 




13 


2 


37 


3 


3 


73 


7 




14 


2 


37 


3 


3 


7 


37 




15 


2 


3 


73 


37 


3 


7 




16 


2 


3 


73 


3 


73 


7 




17 


2 


3 


73 


3 


7 


37 




18 


2 


3 


7 


3 


37 


37 




19 


2 


373 


3 


7 


3 


7 




20 


2 


3 


733 


7 


3 


7 




21 


2 


3 


7 


337 


3 


7 




22 


2 


3 


7 


3 


373 


7 




23 


23 


73 


37 


3 


7 






24 


23 


73 


3 


23 


7 






25 


23 


73 


3 


7 


37 






26 


23 


7 


3 


37 


37 






27 


2 


37 


3 


37 


37 






28 


2 


3 


73 


37 


37 






29 


2 


3733 


7 


3 


7 






30 


2 


3 


7 


3373 


7 






31 


2 


373 


37 


3 


7 






32 


2 


373 


3 


73 


7 






33 


2 


373 


3 


7 


37 






34 


2 


373 


3 


7 


37 






35 


2 


3 


733 


73 


7 






36 


2 


3 


733 


7 


37 






37 


2 


3 


7 


337 


37 






38 


23 


733 


7 


3 


7 






39 


23 


7 


337 


3 


7 






40 


23 


7 


3 


373 


7 






41 


2 


37 


337 


3 


7 
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Table 5. Continued. 



42 


2 


37 


3 




7 


43 


2 


3 


73 


373 


7 


44 


2 


373 


373 


7 




45 


2 


37337 


3 


7 




46 


2 


337 


37 


37 




47 


23 


733 


73 


7 




48 


23 


733 


7 


37 




49 


23 


7 


337 


37 




50 


2 


37 


337 


37 




51 


23 


73 


373 


7 




52 


2 


3733 


73 


7 




53 


2 


3733 


7 


37 




54 


23 


7 


3373 


7 




55 


2 


37 


3373 


7 




56 


2 


37337 


37 






57 


23 


73373 


7 






58 


23 


7 


33737 






59 


2 


37 


33737 






60 


2373373 


7 









With a slight exception for 6-digit integers the ratio between primes that can 
be partitioned into primes and the number of primes is increasing and 
indicates strongly that the prime-partial digital sequences is infinite. 



To list the actual partitions would be impossible but for each n-digit integer 
sequence of primes there are a few primes that are record holders in the 
number of ways in which they can be partitioned into primes (p ways). Here 
is a list of those primes. 



n=3 


P=3 


373 


n=4 


p=5 


3137, 3373, 3733 and 3797 


n=5 


p=12 


37337 


n=6 


p=18 


237373, 537373 and 733373 


n=7 


p=37 


2373373 


n=8 


p=60 


23733737 



We see that of the prime digits 2,3,5 and 7 the primes 3 and 7 play a 
dominant role in the composition of those primes that lend themselves to the 
maximum number of partitions into primes. All partitions can be displayed in 
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explicit form by executing a UBASIC program based on the partition patterns 
which were subject to a detailed discussion in the introduction. Table 5 
shows the partitions for 23733737. 



13 The Saarandache Square-Partial Digital Sequence 

Definition: The square-partial digital sequence is the sequence of perfect 
squares which can be partitioned into two or more groups of digits which are 
also perfect squares. 

The first one hundred terms of this sequence is shown in table 6. 

Table 6. The first one hundred terms f the square-partial cfgital sequence. 



49 


100 


144 


169 


361 


400 


441 


900 


1225 


1369 


1444 


1600 


1681 


1936 


2500 


3249 


3600 


4225 


4900 


6400 


8100 


9025 


9409 


10000 


10404 


11025 


11449 


11664 


12100 


12544 


14161 


14400 


14641 


15625 


16641 


16900 


19044 


19600 


22500 


25600 


28900 


32400 


36100 


36481 


40000 


40401 


41616 


42025 


43681 


44100 


44944 


48400 


49729 


52900 


57600 


62500 


64009 


67600 


72900 


78400 


81225 


84100 


90000 


93025 


93636 


96100 


99225 


102400 


105625 


108900 


115600 


116964 


117649 


119025 


121104 


122500 


129600 


136161 


136900 


140625 


143641 


144400 


152100 


157609 


160000 


161604 


164025 


166464 


168100 


170569 


176400 


184900 


193600 


194481 


202500 


211600 


220900 


225625 


230400 


237169 



This sequence is infinite. If all the infinitely many squares of the form s-lO 2 *, 
where s is a perfect square, are removed the sequence is still infinite. Charles 
Ashbacher [3] has proved this in the case where leading zeros are accepted. 
In order to continue to use unique representation of integers the proof has to 
be changed. 

Proof: There are infinitely many squares of the form (10 ac+1 +2) 2 . In expanded 
form these can be written (10 k ) 2 -10 2<k+1) +(210 k ) 2 -10 +2 2 . From this we see 
that (10 2t,l +2) 2 can be partitioned into the three squares: lO 2 *, -MO 2 * and 4. 
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The methods and programs developed in the introduction to this chapter have 
been used in the analysis of this sequence with the only change that we study 
a different partition property defined though the function 8(t) which now 
takes the form: 

10 .D(t) 

20 if (isqrt(t)) A 2=t then z=l else z=0 
30 retum(z) 

Use has been made of the Ubasic function isqrt(x) which gives the integer 
part of the square root of x. 

The analysis was carried out for 2, 3, ... 8 digit squares. The results are 
presented in the form used in the previous section. 

Table 7. The ratio m/q. m is the number of squares which can be partitioned into 
squares and q is the number of rvcfigit squares. 



n=2 


n=3 


n=4 


n=5 


n=6 


n=7 


n=8 


m 1 


7 


15 


44 


114 


316 


883 


q 6 


15 


68 


217 


683 


2163 


6837 


m/q 0.1667 0 


.4667 0 


.2206 0 


.2028 0 


. 1669 C 


1. 1461 0 


'.1292 




Diagram 5. The ratio m/q for n-digit integers (n - 1,2 ... 8). 
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For squares with more than three digits the graph shows an interesting 
decrease as if it were approaching an horizontal asymptote. Is this the case? 
For each n-digit integer sequence of squares we list those squares which can 
be partitioned into squares in a maximum number of ways (p ways): 



n=3 


p=i 


100, 144, 169, 361, 


400, 441, 900 


n=4 


p=3 


4900 




n=5 


p=5 


11449, 14400, 16900, 


36100, 44100 


n=6 


p=5 


144400, 490000 




n=7 


p=7 


1144900 




n=8 


p=8 


14440000, 36144144, 


49491225 



This list is dominated by squares ending on an even of zoos. The 

digit 4 occurs most frequently while the digits 7 and 8 are lacking. This 
section is terminated by a table of maximum partitions fix n=8. 



Table 8a. Al possibte partitions of 49491 2 25 into sq uares. 



1 


4 


9 


4 


9 


1 


225 


2 


4 


9 


4 


9 


1225 




3 


49 


4 


9 


1 


225 




4 


4 


9 


49 


1 


225 




5 


49 


49 


1 


225 






6 


49 


4 


9 


1225 






7 


4 


9 


49 


1225 






8 


49 


49 


1225 










Table 8b. 


Al possible partitions of 361 441 44 into squares. 




1 


36 


1 


4 


4 1 


4 


4 


2 


361 


4 


4 


1 4 


4 




3 


36 


144 


1 


4 4 






4 


36 


1 


441 


4 4 






5 


36 


1 


4 


4 144 






S 


361 


441 


4 


4 






7 


361 


4 


4 


144 






8 


36 


144 


144 
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Table 8c. All possible partitions of 1 44440000 into squares. 



1 


1 


4 


4 


4 


0 


0 


2 


144 


4 


0 


0 


0 


0 


3 


1 


4 


4 


400 


0 


0 


4 


1444 


0 


0 


0 


0 




5 


144 


400 


0 


0 






6 


1 


4 


4 


40000 






7 


144400 


0 


0 








8 


144 


40000 











L4 The Smarandache Cube-Partial Digital Sequence 

Definition: The cube-partial digital sequence is the sequence of cubes which 
can be partitioned into two or more groups of digits which are also cubes. 

It is obvious that all integers of the form n 3 *10 3k where n is an integer belong 
to the sequence, but if we let n be an integer whose last digit *0, then there 
are only two cubes m=n 3 <10 9 which can be partitioned into cubes. They are: 

22 3 = 10648 with the partition 1, 0, 64, 8, 
and 

303 3 =2781812? with the partition 27, 8, 1, 8, 1, 27. 

In spite of this the sequence has infinitely many members with the last digit 
^0. This is due to the fact that all cubes of the form (3-10 3k+2 +3) 3 can be 
written in the form 27- 10^+8 M0 61c+4 +8M0 3k+2 +27 which gives the 
partition: 

27, 0, 0 ...(3k zeros), 81, 0 ,0 ,..(3k zeros), 81, 0 ,0 ... (3kzeros), 27. 

The integer 303 can be looked at as the generator of a whole family of 
members of the cube-partial digital sequence. 

Question: Is 303 the only generator that guaranties that there are infinitely 
many members of this sequence with non-zero last digits? 
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L5 ParthioB of {1, 2, 3, ... 2b} iato trigrades using Pythagorean triples. 



The relation 2 k +3 k +7 k =l k +5 k +6 k bolds for k=1.2. It is called a bigrade and is 

written 2,3,7 = 1,5,6. To find similar relations which bold for a larger range of 
values for k has attracted a lot of interest [4], Here is an example of a 
pentagrade 

1,5,10,18,23,27 = 2,3,13,15,25,26 

Definition: A relation 

ai m +a2*+ ... bi m +b2 M + ... W* 

which holds for m=I, 2, .... k is called a k-grade. In abbreviated form this is 
written 

^l»32»"*an = bl,b2,...bn 

The purpose of this study is to give an explicit expression for the partition of 

S={1,2, .. 2n} into a={ ai, a^...^} and B={ bj, b 2 ,...,bn} so that A= B \ We 
will see that this partition can be obtained through a suitable choice erf 
Pythagorean numbers. The trigrades generated in this way will in general 
have a large number of terms. They can, however, be used to produce 
trigrades with only eight terms: 

,^ 2,33 a4 = bj,b2,b3,b4 

Let Ao consist of all odd integers x in the interval l^x<n-l, where n is 
assumed to be even, and all even integers in the interval n+2£x^2n. Bo=A 0 
Consider 

S m =— £k m (Half the sum erf m* powers erf S) 

2k=i 



nJ 2 a/2 

A b = I(2k~l) m +I(n+2k) n 

k=l k=l 

Evaluating these sums result in 



1 S = A v-/ B with A n B = 0, b = a’ (a’ stands for the complement of a). 
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( 1 ) 



A, - S, = 0 (i.e. Ao = B 0 ) 

A 2 - S2 = n 2 /2 (2) 

A 3 - S 3 = 3n 2 (2n+l)/4 (3) 

A solution is attempted by constructing new partitions A and B through 
exchange of elements between Ao and Bo so that (1) is maintained while 
satisfying (2) and (3). 

Exchange: (xi , x 2 } e Bo A and {yi , y 2 } e Ao -> B so that 

yi+y 2 = x, +x 2 (1.1) 

(yi 2 + yz 2 ) - (X1 2 + x 2 2 ) = n 2 /2 (2.1) 

(yi 3 + y 2 3 ) - (X , 3 + x 2 3 ) = 3n 2 (2n+l)/4 (3.1) 

Assume y 2 >yi and write (1.1) as: 

xi=yi+d; x 2 = y 2 -d (1.2) 

Substituting this and simplifying (2.1) and (3.1) results in: 

2d(y 2 - yi -d) = n 2 /2 ( 2.2) 

y, +y 2 =2n + l (3 2) 

Assume j s (1, 2, ... n/2} and write (3.2) as 

y 2 = n + 2j; yi = n + 1 - 2j (3.3) 

Insert these expressions in (2.2) 

n 2 = 4d(4j - d -1) (2.3) 
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Integer solutions for d require the discriminant of this equation to be a square 

z 2 . 



(4j - l) 2 - n 2 = z 2 with n/4 + 1 5 n/2 (4) 

Let a, b and c be a Pythagorean triple with 2 lb and let k be an odd integer, 
(4) is then satisfied by 4j-l=kc, n=kb and z=ka. It remains to evaluate xj, x 2 , 
yi and y 2 in terms of a, b, c and k. To begin with we have j=(kc+iy4. Since j 
s Z+ and c=l (mod 4) we can write k=4e-l where e e Z+. This gives d=(4e- 
lXc ± a) and the final solution: 

xu = ((4e-lX2b ± a) + iy2 (5) 

yu = ((4e-lX2b ± c)+iy2 (6) 

The conditions 2b>a and 2b>c restrict the number of Pythagorean triples 
which can be used to generate trigrade partitions of the first N positive 
integers where N is given by: 

N = 2b(4e - 1), e £ Z+ (7) 

A further condition is that xi and x 2 given by (5) must belong to Ao, i.e. xaO 
(mod 2) if CKx^n and x=l (mod 2) if n<x<2n This is met if 2b+a=-l (mod 4). 
The corresponding condition must be met for yi and y 2 . 

From the above conditions and (7) we see that a Pythagorean triple generates 
zero or infintely many trigrades. If it produces infinitely many tri grades then 
the number of terms N of these are in arithmetic progression with first term 
6b and difference 8b. The first non-trivial Pythagorean triples are listed in 
table 9. Those which can be used to produce trigrades are marked in column 
N by giving the number of trams in the trigrade which has the least number of 
terms. 

The trigrade with the smallest number of terms generated in this way is (the 
exchange terms are underscored); 

1, 2, 7, 8, 9, 11, 14, 16, 17, 18, 22, 24 = 2, 4, 5, 6, 10, 12, 13, 15, 19, 20, 21, 23 
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Table 9. Pythagorean triples and their corresponding smallest N. 



a 


b 


c 


N 


a 


b 


c 


N 


a 


b 


c 


N 


3 


4 


5 


24 


55 


48 


73 


288 


17 


144 


145 




5 


12 


13 




77 


36 


85 




51 


140 


149 


840 


15 


8 


17 




13 


84 


85 




85 


132 


157 




7 


24 


25 


144 


39 


80 


89 


480 


119 


120 


169 


720 


21 


20 


29 




65 


72 


97 




165 


52 


173 




35 


12 


37 




99 


20 


101 




19 


180 


181 


1080 


9 


40 


41 




91 


60 


109 


360 


57 


176 


185 




45 


28 


53 




15 


112 


113 


672 


153 


104 


185 




11 


60 


61 


360 


117 


44 


125 




95 


168 


193 


1008 


33 


56 


65 




105 


88 


137 




195 


28 


197 




63 


16 


65 




143 


24 


145 




133 


156 


205 





Since the sets Ao and Bo are defined from the outset the trigrade partitions can 
be unambigously described in terms of these sets and the exchange terms as 
shown in table 10. From this table it is seen that the trigrades corresponding 
to the last two Pythagorean triples have the same number of terms. This 
means that we can obtain trigrades with only eight terms by subtraction. For 
example we see from table 10 that 

197, 164, 344, 17 = 3 17, 44, 272, 89 (8) 

A k-grade a 1 ,a 2 ,...a n = bi,b 2 ,...bn is invariant under translation. To prove this 
consider the sets A={x+a u x+a 2 ,...,x+a n } and B={x+b u x+b^.^x+bn}. 
Expand each term to the power m where m<k 



Z( x + *i) 
i-1 



■zz 



|x J X." 




n * ( m) . m 

sSiK 



i=l 






where the assumed k-grade property has been used in exchanging ^ and bj. 
We can consequently reduce all terms in the above trigrade by 17 (if we 
reduced by 16 instead the eight terms would be positive and realively prime 
but would prevent further reductions) to obtain 



0, 147, 180, 327 = 27, 72, 255, 300 



( 8 . 1 ) 
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The k -grade property is evidently invariant under division by a common 
factor. Hence 

0, 49, 60, 109 = 9, 24, 85, 100 (8.2) 

The steps shown above carried out for Pythagorean triples c<1000 result in 
quite a large number of eight-term trigrades as is shown in table 11. To 
reduce the trigrades so that the smallest term is 0 and divided the terms with 
their LGM proved essential in order to avoid duplications. When 
implementing the algorithms on a computer several solutions woe obtained 
which after reduction proved to be identical. 

Table 10. Description of trigrades 



# 

T 

2 

3 

4 

5 

6 

7 

8 
9 



Tripl« 

V 



a-*B 



3,4,5 

24 

56 

88 

7,24,25 

144 

336 

528 

55, 48,73 
288 

672 

1056 



a~{l,3,-,li}o{14,16,_,24}w{8,17}-{5 f 20 } 

B-{2,4,-,12}u{13,l5,_,23}u{5,20}-{8,17} 

A~{1,3,_,27}vj{30,32,~,56}w{18,39)-{11,46} 

B«{2,4,_,28}o{29,31,_,55}w{U,46}-{18,39} 

A~{1,3,-,43}w{46,48,_,88}w{28,61}-{17,72} 

B-{2,4,_,44)w{45,47,_,87}^{17,72}-{28,61} 

A={1,3,_,71}w{ 74,76,„, 144}w{62,83}-{35,110} 

B={2,4,_,72}w{73,75,_,143}w{35,110}-{62,83) 

AH1,3 # -.,167}w{170,172,~,336}w{144, 193} -{ 81 , 256} 

B={2,4,_,168}w{169,171,„,335}^{81,256}-{144,193} 

A={1,3,-.,263}w{266,268,„,528}^{226, 303} -{ 127, 402 } 

B={ 2, 4, ~, 264 } w{ 265, 267 527 }w{ 127, 4 02 }-{ 22 6, 303} 

A=*{l,3,_,l43}w{146,148,~,288}u>{62,227}-{35,254} 

B-{2,4,-,144}o{145,147,„,287}o{35,254}-{62,227} 

A={ 1, 3,~, 335}v{338, 34 0 ,-., 672 } 144 , 529} -{ 81, 592 } 

B={ 2, 4,„, 336} < j { 337, 339, 671 }w{ 81, 592 >-{144,529} 
A={l,3,_,527}w{530,532,_,1056}^{226,831}-{127, 930} 
B={2,4,„ r 528}w{52 9,531,^.,1055}w{127, 930} -{ 226, 831 } 



Table 10. Description of trigrades, continued. 



♦ Tr± * lm a^B 

H 

10 11,60,61 A={l,3,„,179}w{182,184,~,360}u{164,197}-{89,272} 

360 B={ 2, 4,_, 180} w{ 181, 183, 359} w{ 8 9, 272} -{164, 197} 
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11 




A= ( 1, 3, 41 9)^(422, 424, -, 




840 


B= ( 2, 4,_,420}u{ 421, 4 23,-, 


12 


91,60, 


109 A=( 1 ,3,— ,179)vj(182,184,.... 




360 


B={2,4,_,180}w(181,183,_, 


12 




A=(1,3,._,419)vj{422,424,..., 




840 


B=( 2, 4,—, 420} V— » { 421 , 423,—, 



840)^(382,459} -{207, 634} 

839}w{207,634}-{382,459} 

360}w{44,317}-{17,344} 

359<~<{17,344)-(44,317} 

840}w{102,739}-(39,802} 

839}^(39, 802} -(102, 739} 



Table 1 1 . Trigrades 



Si 


a 2 


a 3 


a 4 


bi 


b 2 


b 3 


b 4 


0 


49 


60 


109 


9 


24 


85 


100 


0 


67 


102 


169 


22 


25 


144 


147 


0 


85 


136 


221 


25 


36 


185 


196 


0 


103 


122 


225 


22 


45 


180 


203 


0 


107 


242 


349 


25 


62 


287 


324 


0 


125 


638 


763 


50 


63 


700 


713 


0 


129 


244 


373 


48 


49 


324 


325 


0 


151 


382 


533 


49 


74 


459 


484 


0 


155 


230 


385 


34 


77 


308 


351 


0 


165 


184 


349 


25 


84 


265 


324 


0 


181 


384 


565 


60 


81 


484 


505 


0 


185 


262 


447 


10 


147 


300 


437 


0 


203 


554 


757 


81 


86 


671 


676 


0 


215 


370 


585 


46 


117 


468 


539 


0 


229 


304 


533 


49 


108 


425 


484 


0 


233 


278 


511 


63 


86 


425 


448 


0 


241 


556 


797 


72 


121 


676 


725 


0 


271 


678 


949 


102 


117 


832 


847 


0 


283 


542 


825 


58 


165 


660 


767 


0 


293 


370 


663 


75 


118 


545 


588 


0 


301 


456 


757 


81 


132 


625 


676 


0 


305 


458 


763 


63 


158 


605 


700 


0 


331 


994 


1325 


110 


169 


1156 


1215 


0 


343 


1074 


1417 


117 


174 


1243 


1300 






Table 1 1 


, Trigrades, continued. 








a l 


a 2 


a 3 


a 4 


bi 


b 2 


b 3 


b 4 


0 


349 


376 


725 


49 


180 


545 


676 


0 


359 


506 


865 


81 


170 


695 


784 


0 


365 


400 


765 


76 


153 


612 


689 


0 


379 


522 


901 


54 


225 


676 


847 


0 


381 


640 


1021 


121 


156 


865 


900 



25 








0 


385 


996 


1381 


96 


225 


1156 


1285 


0 


407 


1262 


1669 


122 


225 


1444 


1547 


0 


461 


1392 


1853 


153 


236 


1617 


1700 


0 


469 


1264 


1733 


108 


289 


1444 


1625 


0 


477 


764 


1241 


153 


188 


1053 


1088 


0 


501 


688 


1189 


76 


289 


900 


1113 


0 


511 


1348 


1859 


148 


275 


1584 


1711 


0 


541 


736 


1277 


121 


252 


1025 


1156 


0 


551 


796 


1347 


147 


236 


1111 


1200 


0 


593 


824 


1417 


121 


296 


1121 


1296 


0 


631 


1138 


1769 


169 


306 


1463 


1600 


0 


649 


964 


1613 


169 


288 


1325 


1444 


0 


787 


838 


1625 


162 


325 


1300 


1463 


0 


817 


1048 


1865 


245 


292 


1573 


1620 


0 


835 


1880 


2715 


296 


363 


2352 


2419 


0 


845 


2522 


3367 


170 


567 


2800 


3197 


0 


849 


1084 


1933 


169 


408 


1525 


1764 


0 


925 


2190 


3115 


34 


847 


2266 


3081 



L6 Smarandache Non-Null Squares. 

Question: In how many ways can n be written as a sum of non-null squares, 
disregarding the order of the terms? 

Example: 9=l 2 +l 2 +l 2 +l 2 +l 2 +l 2 +l 2 +l z +l 2 

= 1 2 + 1 2 + 1 2 + 1 z + 1 2 + 2 2 
= 1 2 + 2 2 + 2 z 
=3 2 

Let us denote the number of ways by i?(n). The example shows that 
r)(9>=4. Obviously the non-null squares representation which has the largest 
number of terms is a r ep re se n tation by n squares of 1 ami the representation 
with the smalle st number of terms is a representation where we use the largest 
possible squares. If n is a perfect square then (Vn ) 2 is one -term 
representation. 

Let’s assume that we know rj(x) for all integers x<n-l. Furthermore we make 
the definition: 

. . J 1 if n is a square 
(0 otherwise 



26 





With these precisions we can write the algorithm for calculation of r|(n) in 
the following way: 



t)(n) = l+8(n) 



[Vn]-5(n) 

+ z 

1=2 






with r|(l) = 1 



This algorithm was implemented for n<100. As is shown in table 12 this 
sequence, after a slow start, grows very rapidly. 



Table 12. The Non-Nufl Square Representations. 



1 


1 


1 


2 


2 


2 


2 


3 


4 


4 


4 


5 


7 


7 


7 


9 


12 


13 


13 


16 


20 


23 


23 


27 


35 


41 


42 


47 


61 


71 


75 


82 


104 


124 


134 


146 


178 


217 


237 


258 


307 


377 


419 


456 


535 


651 


739 


804 


933 


1126 


1300 


1422 


1629 


1955 


2275 


2513 


2846 


3397 


3972 


4435 


4990 


5904 


6933 


7 807. 


8766 


10268 


12097 


13718 


15409 


17895 


21087 


24076 


27076 


31248 


36736 


42214 


47568 


54636 


64017 


73924 


83554 


95596 


111637 


129306 


146714 


167379 


194807 


226021 


257447 


293255 


340106 


394953 


451408 


514025 


594103 


690035 


790967 


901118 


1038450 


1205451 



L7 Smarandache Non-Null Cubes. 



In how many ways x(n) can n be written as a sum of non-null cubes? In this 
case there are much fewer representations. For non-null squares we have 
r|(9)=4 whereas for non-null cubes we have t( 9)=2. Calculation of the first 
100 terms follows the same method as for the non-null squares with 
corresponding changes in our algorithms. 



5{n) = j lif 
^ 1 [Oot 



nisa cube 
otherwise 



x(n) = l+8(n) 



+ 



[l/n]-5(n) 

i 




with x(l) = 1 



Table 13. The Non-NuU Cube Representations. 
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1 


1 


1 


1 


1 


1 


1 


2 


2 


2 


2 


2 


2 


2 


2 


3 


3 


3 


3 


3 


3 


3 


3 


4 


4 


4 


5 


5 


5 


5 


5 


6 


6 


6 


8 


8 


8 


8 


8 


9 


9 


9 


12 


12 


12 


12 


12 


13 


13 


13 


17 


17 


17 


18 


18 


19 


19 


19 


24 


24 


24 


27 


27 


29 


29 


29 


35 


35 


35 


41 


41 


44 


44 


44 


51 


51 


51 


61 


61 


65 


66 


66 


74 


74 


74 


89 


89 


94 


98 


98 


109 


109 


109 


130 


130 


136 


146 


146 


162 


162 
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Chapter II 

Recursive Integer Sequences 



IL 1 The Non-Arithmetic Progression 

We consider an ascending sequence of positive integers ai, a 2 , ... an such that 
each element is as small as possible and no t-term arithmetic progression is in 
the sequence. In order to attack the problem of building such sequences we 
need a more operational definition. 

Definition: The t-term non-arithmetic progression is defined as the set : 

{a, :ai is the smallest integer such that a 1 >a l . 1 and such that for k<i there are at 
most t-1 equal differences a. — a, = a, - a, = . - a, -a. \ 

K *1 *1 *2 2 k t - 1 } 

From this definition we can easily formulate the starting set of a t-term non- 
arithmetic progression: 

{h 2, 3 t-1, t+1} or {at : af=i for i<t-l and a*=t+l where t^3} 

It may seem clumsy to bother to express these simple definitions in stringent 
terms but it is in fact necessary in order to formulate a computer algorithm to 
generate the terms of these sequences. 

Question: How does the density of a t-term non-arithmetic progression vary 
with t, i.e. how does the fraction a^/k behave for t>3? 

Strategy for building a t-term non-arithmetic progression: Given the 

terms ai, a 2 , ... at we will examine in turn the following candidates for the 
term a^+ 1 : 
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aw = a* + d, d=l, 2, 3, ... 



Our solution is the smallest d for which none of the sets 



{ai, a 2 , ... at, at+d, ak+d-e, %+d-2e, ... ak+d-(t*l)-c :esd} 



contains a t-tenn arithmetic progression 

We are certain that a^j exists because in the worst case we may have to 
continue constructing sets until the term ak+d-(t-l)e is less than 1 in which 
case all possibilities have been tried with no t terms in arithmetic progression. 
The method is illustrated with an example in diagram 1. 

In the computer application erf the above method the known terms of a no t- 
term arithmetic pr ogres si on were stored in an array. The trial terms were in 
each case added to this array. In the example we have for d=l, e=l the array: 
1,2,3,5,6,8,9,10,11,10,9,8. The terms are arranged in ascending order 
1,2,3,5,6,8,8,9,9,10,10,11. Three terms 8,9 and 10 are duplicated and 11 
therefore has to be rejected. For d=3, e=3 we have 1,2,3,5,6,8,9,10,13,10,7,4 
or in ascending order 1,2,3,4,5,6,7,8,9,10,10,13 this is acceptable but we 
have to check for all values of e that produce terms which may form a 4-term 
arithmetic progression and as we can see from diagram 1 this happens for 
d=3, e=4, so 13 has to be rejected. However, for d=5, e=5 no 4-tom 
arithmetic progression is formed and e=6 does not produce terms that need to 
be checked, hence a? = 15. 







1 


2 3 


4 5 


6 


7 8 


9 


10 11 


12 


13 


14 


15 




Known terms 




1 


2 3 


5 


6 


8 


9 


10 












Triab 




























d=l 


e=t 










8 


9 


10 11 










reject 1 1 


d=2 


e=2 








6 


8 




10 


12 








reject 12 


d=3 


e=3 






4 




7 




10 




13 






try next e 




e~4 


1 




5 






9 






13 






reject 13 


d=4 


e*3 




2 




6 






10 






14 




reject 1 4 


d=5 


e=5 






5 








10 








15 


accept 15 



Diagram 1 . To find the 9* term of the 4-term non-arithmetic progression. 
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Routines for ordering an array in ascending order and checking for 
duplication of terms were included in a QBASIC program to implement the 
above strategy. 




Diagram 2. a*/k for non-arithmetic progressions with t=3, 4, 5, ... 1 5. Bars are shown 
for k = multiples of 10. 

Results and observations: Calculations were carried out for 3<t<15 to find 
the first 100 terms of each sequence. The first 65 terms and the 100** term are 
shown in table 1. In diagram 2 the fractions a^/k has been chosen as a 
measure of the density of these sequences. The looser the terms are packed the 
larger is ajk. In fact for t>100 the value of ak/k = 1 for the first 100 terms. 




In table 1 there is an interesting leap for t=3 between the 64“ and the 65“ 
terms in that = 365 and - 730. Looking a little closer at such leaps we 
find that: 



Table 1 . The 65 first terms of the norvarithmefic progressions for t=3 to 15. 



# 


♦-3 


t=4 


♦=5 


♦=6 


t*7 


t=8 


t=9 


t=10 


MI 


t*12 


♦=13 


♦=14 


♦=15 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


3 


4 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


4 


5 


5 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


5 


10 


6 


6 


5 


5 


5 


5 


5 


5 


5 


5 


5 


5 


6 


11 


8 


7 


7 


6 


6 


6 


6 


6 


6 


6 


6 


6 


7 


13 


9 


8 


8 


8 


7 


7 


7 


7 


7 


7 


7 


7 


8 


14 


10 


9 


9 


9 


9 


8 


8 


8 


8 


8 


8 


8 


9 


28 


15 


11 


10 


10 


10 


10 


9 


9 


9 


9 


9 


9 


10 


29 


16 


12 


12 


11 


11 


11 


11 


10 


10 


10 


10 


10 


11 


31 


17 


13 


13 


12 


12 


12 


12 


12 


11 


11 


11 


11 


12 


32 


19 


14 


14 


13 


13 


13 


13 


13 


13 


12 


12 


12 


13 


37 


26 


16 


15 


15 


14 


14 


14 


14 


14 


14 


13 


13 


14 


38 


27 


17 


17 


16 


16 


15 


15 


15 


15 


15 


15 


14 


15 


40 


29 


18 


18 


17 


17 


16 


16 


16 


16 


16 


16 


16 


16 


41 


30 


19 


19 


18 


18 


17 


17 


17 


17 


17 


17 


17 


17 


82 


31 


26 


20 


19 


19 


19 


18 


18 


18 


18 


18 


18 


18 


83 


34 


27 


22 


20 


20 


20 


20 


19 


19 


19 


19 


19 


19 


85 


37 


28 


23 


22 


21 


21 


21 


20 


20 


20 


20 


20 


20 


86 


49 


29 


24 


23 


23 


22 


22 


21 


21 


21 


21 


21 


21 


91 


50 


31 


25 


24 


24 


23 


23 


23 


22 


22 


22 


22 


22 


92 


51 


32 


26 


25 


25 


24 


24 


24 


24 


23 


23 


23 


23 


94 


53 


33 


33 


26 


26 


27 


25 


25 


25 


24 


24 


24 


24 


95 


54 


34 


34 


27 


27 


28 


26 


26 


26 


25 


25 


25 


25 


109 


56 


36 


35 


29 


28 


29 


27 


27 


27 


27 


26 


26 


26 


no 


57 


37 


36 


30 


30 


30 


28 


28 


28 


28 


28 


27 


27 


112 


58 


38 


37 


31 


31 


31 


31 


29 


29 


29 


29 


28 


28 


113 


63 


39 


39 


32 


32 


32 


32 


X 


30 


30 


30 


29 


29 


118 


65 


41 


43 


33 


33 


33 


33 


31 


31 


31 


31 


31 


30 


119 


66 


42 


44 


34 


34 


34 


34 


32 


32 


32 


32 


32 


31 


121 


67 


43 


45 


36 


35 


37 


35 


34 


33 


33 


33 


33 


32 


122 


80 


44 


46 


37 


37 


38 


36 


35 


35 


34 


34 


34 


33 


244 


87 


51 


47 


38 


38 


39 


37 


36 


36 


35 


35 


35 


34 


245 


88 


52 


49 


39 


39 


40 


38 


37 


37 


36 


36 


36 


35 


247 


89 


53 


50 


40 


40 


41 


39 


38 


38 


37 


37 


37 


36 


248 


91 


54 


51 


41 


41 


43 


41 


39 


39 


38 


38 


38 


37 


253 


94 


56 


52 


50 


42 


44 


42 


40 


40 


40 


39 


39 



Table continued on next page. 
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Table 1 continued 




Leap starts at Leap finishes at 



5 




10 




14 


= 3 - 5-1 


28 


= 2-14 


41 


= 3 - 14-1 


82 


= 2-41 


122 


= 3 - 41-1 


244 


= 2-122 


365 


= 3 - 122-1 


730 


= 2-365 



Does this chain of regularity continue indefinitely? 
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Sometimes it is easier to look at what is missing than to look at what we have. 
Here are some observations on the only excluded integers when fo rming the 
first 100 terms for t= 11, 12, 13 and 14. 

For t=ll: 11, 22, 33, 44, 55, 66, 77, 88, 99 The nth missing integer 

is ll*n 

For t=12: 12, 23, 34, 45, 56, 67, 78, 89, 100 The n th missing integer 

is 11-n+l 

For t=13: 13, 26, 39, 52, 65, 78, 91, 104 The n th mi ssing integer 

is 13-n 



For t=14: 14, 27, 40, 53, 66, 79, 92, 105 The nth missing integer 

is 13-n+l 



Do these regularities erf missing integers continue indefinitely? What about 
similar observations for other values erf t? 



n.2 The Prime-Product Sequence 

Definition: The terms erf the prime-product sequence are defined through {to : 
to = p»#+l, Pn is the nth prime number}, where Po# denotes the product of all 
prime numbers which are less than or equal to Pn. 

The sequence begins {3, 7, 31, 211, 2311, 30031, ... }. In the initial 
definition of this sequence ti was defined to be equal to 2. However, there 
seems to be no reason for this exception. 

Question: How many members erf this sequence are prime numbers? 

The question is in the same category as questions like How many prime twins 
are there?. Haw many Carmichael numbers are there?, etc. ’ So we may have 
to contend ourselves by finding how frequently we find prime numbers when 
examining a fairly large number of terms of this sequence. 
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From the definition it is clear that the smallest prime number which divides tn 
is larger than pn. The terms of this sequence grow rapidly. The prime number 
functions prmdiv(n) and nxtprmfn) built into the Ubasic programming 
language were used to construct a prime factorization program for n<10 19 . 
This program was used to factorize the 18 first terms of the sequence. An 
elliptic curve factorization program, ECM.UB, conceived by Y. Kida was 
adapted to generate and factorize further terms up to and including the 49 th 
term. The result is shown in table 2. All terms analysed were found to be 
square free. A scatter diagram. Diagram 3, illustrates how many prime factors 
there are in each term . 

The 50 th term presented a problem. t 5 o= 126173 *n, where n has at least two 
factors. At this point prime factorization begins to be too time consuming and 
after a few more terms the numbers will be too large to handle with the above 
mentioned program. To obtain more information the method of factorizing 
was given up in favor of using Fermat’s theorem to eliminate terms which are 
definitely not prime numbers. We recall Fermat’s little theorem: 

If p is a prime number and (a, p)=l then a**" 1 * 1 (mod p). 

a 0-1 = 1 (mod n) is therefore a necessary but not sufficient condition for n to be 
a prime number. If n fills the congruence without being a prime number then 
n is called a pseudo prime to the base a, psp(a). We will proceed to find all 
terms in the sequence which fill the congruence 

a tm ~ 1 = l(modf rt ) 

for 50 <,n<, 200. t 2 oo is a 5 13 digit number so we need to reduce the powers of 
a to the modulus tn gradually as we go along. For this purpose we write tn-1 to 
the base 2: 



to-l = S(k)'2 k , where 5(k) e {0,1} 

k=\ 

From this we have 

Jr=l 
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Table 2* Prime factorization of prime-product terms 



a 


D 


D 


N=p#+1 and rh focton 


i 


2 


1 


3 Prim* number 


2 


3 


1 


7 Prime number 


3 


5 


2 


31 Prime number 


4 


7 


3 


21 1 Prime number 


5 


11 


4 


231 1 Prime number 


6 


13 


5 


30031 =59-509 


7 


17 


6 


510511 =19-97-277 


8 


19 


7 


9699691 =347 * 27953 


9 


23 


9 


223092871 =317 - 703763 


10 


29 


10 


6469693231 =331 - 571 - 34231 


11 


31 


12 


200560490131 Prime number 


12 


37 


13 


7420738134811 =181 - 60611 • 676421 


13 


41 


15 


30425026352721 1 =61 ■ 450451 - 1 1072701 


14 


43 


17 


13082761331670031= 167-78339888213593 


15 


47 


18 


61488978258849141 1 = 953 • 46727 - 13806181 181 


16 


53 


20 


32589158477190044731 = 73 - 139 - 173 - 18564761860301 


17 


59 


22 


1922760350154212639071 =277 - 3467 - 105229 - 19026377261 


18 


61 


24 


1 17288381359406970983271 =223 - 525956867062542470777 


19 


67 


25 


7858321551080267055879091 =54730729297- 143581524529603 


20 


71 


27 


557940830126698960967415391 = 1063 - 303049 - 598841 - 2892214489673 


21 


73 


29 


40729680599249024 1 5062 1 323471 =2521 - 16156160491570418147806951 


22 


79 


31 


3217644767340672907899084554131 =22093- 1503181961 - 96888414202798247 


23 


83 


33 


267064515689275851355624017992791 = 265739 - 1004988035964897329167431269 


24 


89 


35 


2376874189634555077065053760135831 1 = 131 • 1039 - 2719 - 
6422589 1 88429437337 1 806 1 4 1 


25 


97 


37 


2305567963945518424753102147331756071 =2336993- 13848803- 
71 237436024091007473549 


26 


101 


39 


232862364358497360900063316880507363071 = 960703 • 
2423874645530380990795941 27X1 057 


27 




41 


239848235289252281 72706521 63869225839621 1 = 2297 - 9700398839 
1 79365737X7 - 6X1 31 544333453 1 


28 




43 


256637611759499941447959781534X71648394471 = 149 -13203797 - 
30501264491063137 -4276784365106371 1 


29 




45 


279734996817854936178276161872067809674997231 = 334507-1290433- 
6480464442342997146231 77554034701 


30 




47 


31 61 005464041 76077881 4520629 1 54366249327468699 1=51 22427 - 
2025436786X7-30467075950695402471 570558 1 9 


31 




49 


4014476939333036189094441199026045136645885247731 = 
1543-49999-552X1 -57900988201093-1628080529999073967231 
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51 


525896479052627740771371797072411912900610967452631 = 
1951 -22993-1172323185947X14144932345466415143728266617 
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Tabte Z Prime factorization of prime-product terms, continued. 
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Diogrcm 3. The number of prime factors in the first 49 terms of the prime-product 
sequence. 



This product expression for a 1 9 1 is used in the following Ubasic program to 

carry out the reduction erf a^ 1 modulus U, Terms for which 8(k)=0 are 
ignored in the expansion were the exponents k are contained in the array 
E%0 The residue modulus k is stored in F. In the program below the 
reduction is done to base A=7. 

100 dm E%{1000) 

110 M=N-1:I%=0 
120 T=1J%=0 
130 whfle (M-T)>=0 
140 inc J%:T=2*T 
150 wend 

160 dec J%:M=M-T\2:inc l%:E%fl%)=J% 

1 70 if M>0 then goto 120 
180 F=1 

190 for J%=1 tol% 

200 A=7 

210 for K%=] to 

240 A=(AA2)@N 

250 next 

260 F=F*A:F=F@N 

270 next 
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N=20632S54441 86971 805262785592040287445726865285688900747340490078401 81 45 
718728624430191 58728631 60685721 48631 38937930928474301 6940885980671 8870830 
26597753881 31 7772605885038331 62528205231 1 1 21 3067921 93540483321 70364563007 
1 7761 688853571 2671 50232508655634427663661 80331 20098071 1 247645589424056809 
05346832390674579572622346848343362525900088741 1 9591 9732397361 34883450319 
1 30587753586846905761 460662768750585961 002361 1 2260054944287636531 

The last two primes or pseudo primes are remarkable in that they are 
generated by the prime twins 1019 and 1021. 

Summary of results: The number of {nimes q among the first 200 terms of 
the prime -product sequence is given by 6^q<9. The six confirmed primes are 
terms numero 1, 2, 3, 4, 5 and 11. The three terms which are either {Himes or 
pseudo primes are terms numero 75, 171 and 172. The latter two are the 
terms 1019#+1 and 1021#+1. 



1L3 The Square-Product Sequence 



Definition: The toms of the square-product sequence are defined through {tn 
: tn = (n!) 2 +l} 
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This sequence has a structure which is similar to the prime-product sequence. 
The analysis is therefore carried out almost identically to the one done for the 
prime-product sequence. We merely have to state the results and compare 
them. 

The sequence begins {2, 5, 37, 577, 14401, 518401, ... }.A s for the prime* 
product sequ en ce the question of how many are prime numbers has been 
raised and we may never know. There are similarities between these two 
sequences. There are quite a few primes among the first terms. After that they 
become more and more rare. Complete factorization of the 37 first terms of 
the square-factorial sequence was obtained and has been used in diagram 4 
which should be compared with the corresponding diagram 3 for the prime- 
product sequence. 




Diagram 4. The number of prime factors in the first 40 terms of the square-factorial 
sequence. 

Diagram 4 is based cm table 3 which shows the prime factorization of the 40 
first terms in the square-factorial sequence. The number of factors of each 
term is denoted f. The factorization is not complete for terms numero 38 and 
39. A +-sign in the column for f in table 3 indicates that the last factor is not a 
prime. The terms of this sequence are in general much more time consuming 
to factorize than those of die prime-product sequence which accounts for the 
more limited results in this section. Using the same method as for the prime- 
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Table 3. Prime fadorizafion of square-factorial terms. 



N^nlp+l and Hi factors 



l 37 
1 577 

1 14401 

2 518401=13*39877 

2 25401601*101-251501 
2 1625702401= 17- 95629553 
1 131681894401 
1 13168189440001 

1 1593350922240001 

2 229442532802560001=101- 2271708245569901 

1 38775788043632640001 

3 7600054456551997440001* 29- 109- 2404319663572286441 

2 1710012252724199424000001=1344169-1272170577304043929 

2 437763136697395052544000001=149-2938007628841577533852349 
2 126513546505547170185216000001=9049 13980942259426143240713449 
2 40990389067797283140009984000001=37-1107848353183710355135404972973 

2 14797530453474819213543604224000001=710341-20831587158104092560535861261 
5 5919012181389927685417441689600000001=41-10657-86816017* 

348046955609-448324749841 

3 2610284371992958109269091785113600000001= 

6M 57-2725576247251 705240961 774861 76631 51 3 

4 1263377636044591724886240423994982400000001=337 8017-514049836440277481 
909674823323537849 

3 668326769467589022464821184293345689600000001=509-15448374629- 
84994002604532747687401 741 72344 1 
1 384956219213331276939737002152967117209600000001 



941-81 5769831 908479758733- 31 342S331 349331 2902433994 1 7 
5 162644002617632464507038883409628607021056000000000001 = 

53- 53- 418633- 6017159668589- 22985889712876096222556462301797 



113-42461- 745837- 2460281- 7566641- 15238649- U 6793504008451 126962009 
92956902680071 18823449497526840549554238626201 6000000X8001 = 
2122590346576634509- 43794085292997939303952241474982753464389 
78176755153939869305210274200729021751 1 46846355456000000000001 = 
171707860473207588349837- 455289320701414063716469396531 758248773 



I . »’v.l 1 vj M fj Lvl [i' r i mimimmi *TF 



61- 1733- 15661- 359525849- 100636381 126568690110069- 1 174592249518207759537897 



I f*/i w/«'7-v^rA f/.Vj [iVx • 



353-422041- 13400767181- 33867608180948409085305820793832191570324667821677 
6923783734542601 51931661 9894395958365446851 31 90741907865600000000000001 = 
10591621681-6415450638021-522303293914660001204969-1950882388585355532025429 



37- 3121 4421- 4073332882845936253- 3625813512324448042738745076 2106578 223148 
052301 
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Tobte 3 continued. 



□ 


D 


B 


N=M)2+1 and Hi facto n 


34 

35 


77 

81 


4 

3 


871 624056287592837104338371 91 36761 6509271 93975961 39359414779904000000000000 
01=193- 3217- 8661 00731 693394521 434431 4564523 1 476894644096901 291 1 9741 0624286 
816576197 

1 067739468952301 22545281 4505594253302238581 2620552707152831 0538240000000000 
000001=317-3739030196538868064889782854556323553608634748201 176163219890777 
16189815715361 


36 


84 


3 


138379035176218238818684759925015227970120131562363064700690457559040000000 
000000001*73- 57986941373 *3269017431698277804505207628624951 781 7882480631366 
14675401 1637734149869 


37| 


87 


3 


1 89440899 1 5624276894277943633734584709 1 094460 1 08875062955245236398325760000 
000000000001=127406364297881-49105571 194338128021910109- 302797201 1452403842 
829243076981 4272052327643069 


38 


90 


4* 


273552658381614558353373506071 127403199540400397215590907374121359182397440 
000000000000001=233-757-15509190807491 428121 70094885906800637254241 6722731 79 
032363883419184506253167858216021 


39 


93 


4+ 


416073593398435743255481 1027341847802665009490041649137701 16038587316426506 
240000000000000001 *61-1 004545757741 -679001 282693270665470645882 1 250934897329 
61 536897293728551 0229243) 5727303361 7801 




96 


4 


6657177494374971 89208769764374695648426401 51 8406663862032185661 739706282409 
984000000000000000001=89- 701 187100101949-5703061928798691 51956315673142222 
36213934965395443794926477317483490202551 13441 



product sequence the terms t* in the interval 40< n<2 00 which may possible be 
primes were identified- There arc only two of them: 

Tetm »65. N-(651) ? -*-l . N is a 182 doit number. 

68023740289078328950450781 972622203792902576953271 3580342793801 040271 0065 
24643826496596237244465781 5141 2858996571 534385340563792951822384455180747 
800576000000000000000000000000000001 

Term #76 N=(76l) 2 +1 . N is o 223 digit number . 

355509027001 07478542025131 357707726481 94325666925541 647977005250280050084 
1 772266884421 391 665890651 64392091 2930369944999452531 006264950776782697850 
71 9865801 1 625298409931 764786386381 15061 7600000000000000000000000000000000 
0001 



IL4 The Snuurandache Prime-Digital Sub-Sequence 

Definition: The prime-digital sub-sequence is the set : 
{M=ao+ai-10+a 2 -10 2 +...akl0 k Mis prime, all digits ao, ai, a 2 ...% are primes} 
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The first terms of this sequence are {2, 3, 5, 7, 23, 37, 53, 73, . .. >. Sylvester 
Smith [1] conjectured that this sequence is infinite. In this paper we will 
prove that this sequence is in fact infinite. Let's first calculate some more 
terms of the sequence and at the same time find how many terms there is in 
die sequence in a given interval, say between 10* and lO** 1 The program 
below is written in Ubasic . One version of the program has been used to 
produce table 4 showing the first 100 terms of the sequence. The output of the 
actual version has been used to produce the calculated part of table 5 which 
we are going to compare with the theoretically estimated part in the same 
table. 



Ubasic program 



‘Digits not avowed stored in A%() 

'Digits aflowed stored in B%() 
’Calc, for 7 separate intervals 



10 point 2 
20 dferi A%{6),B%(4} 

30 for l%=! to 6:read A%[!%)mext 
40 data 1 ,4,6,8,9,0 
50 for I%=1 to 4read B%(l%):next 
60 dataZ3,5,7 
70 for K%=] to 7 
00 M%=O:N=0 

90 for E%=1 to 4 Only 2. 3, 5, 7 a»owed as first digit 

100 P=B%fE%J*10A|(%:P0=P^=(B%|E%)+1)*10AK%:gosub 150 
110 next 

120 print K%>A%,N,M%/N 
130 next 
140 end 
150 while P<$ 

160 P=nxtpfm(P):P$=str(P} 

170 inc N 

180 L%=)en(P$)£%=0 
190 for f%=2 to L% 

200 for J%=1 to 6 

210 if val(mid(P$,l%,l))=A%{J%) then C%=1 
220 nextr>ext 
230 if C%=0 then inc M% 

240 wend 
250 return 



'Select prime and convert to string 

‘Count number of primes 

l C% wUI be set to 1 if P not member 

‘This loop examines each digit of P 



‘If criteria f*ed count member (m%) 
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Table 4. The first 100 terms in the prime-<Sgital sub sequence. 



2 


3 


5 


7 


23 


37 


53 


73 


223 


227 


233 


257 


277 


337 


353 


373 


523 


557 


577 


727 


733 


757 


773 


2237 


2273 


2333 


2357 


2377 


2557 


2753 


2777 


3253 


3257 


3323 


3373 


3527 


3533 


3557 


3727 


3733 


5227 


5233 


5237 


5273 


5323 


5333 


5527 


5557 


5573 


5737 


7237 


7253 


7333 


7523 


7537 


7573 


7577 


7723 


7727 


7753 


7757 


22273 


22277 


22573 


22727 


22777 


23227 


23327 


23333 


23357 


23537 


23557 


23753 


23773 


25237 


25253 


25357 


25373 


25523 


25537 


25577 


25733 


27253 


27277 


27337 


27527 


27733 


27737 


27773 


32233 


32237 


32257 


32323 


32327 


32353 


32377 


32533 


32537 


32573 


33223 



Table 5. Comparison of results. 



— s — 1 


1 


2 


3 


4 


5 


6 


7 


Computer cotnt 














m 


4 


15 


38 


128 


389 


1325 


4643 


k>g(m} 


0.6021 


1.1761 


15798 


2.1072 


25899 


3.1222 


3.6668 


n 


13 


64 


472 


3771 


30848 


261682 


2275350 


m/n 


0.30769 


0.23438 


0.08051 


0.03394 


0.01261 


0.00506 


0.00204 




■nates: 














m 


4 


11 


34 


109 


364 


1253 


4395 


k>g(m) 


05922 


1.0430 






25615 


3.0980 


3.6430 


n 


7 


55 


421 


3399 


28464 


244745 


2146190 


m/n 

















Theorem: 

The Smarandache prime-digital sub sequence is infinite. 

Proof; 

We recall the prime counting function Jt(x). The number of primes p£x is 
de note 7t(x). For sufficiently large values of x the order of magnitude of 









X 



7i(x) is given by x(x) * 



logx 



Let a and b be digits such that a>b*0 and 



n(a,bjc) the approximate number of primes in the interval (b-lO^a-lO*). 
Applying the {Mime number counting theorem we then have: 



10 * 

n(a,b,k)*—( 



log a logi 

logl0 +— — logl0 + — r— 



) 



( 1 ) 



Potential candidates for members of the prime-digital sub sequence will have 
first digits 2,3,5 or 7, i.e. for a given k they will be found in the intervals 
(2* 10^4- 10*), (5-10^6*10^) and (T-IO^-IO*). The approximate number of 
primes n(k) in the interval (lO^lO 1 ^ 1 ) which might be members of the 
sequence is therefore : 

n(k)=n(4,24cKn(6,5,k)+n(8,7,k) (2) 



The theoretical estimates of n in table 5 are calculated using (2) ignoring the 
fact that results may not be all that good for small values of k. 

We will now find an estimate for the number of candidates m(k) which 
qualify as members of the sequence. The final digit of a prime number >5 can 
only be 1,3,7 or 9. Assuming that these will occur with equal probability only 
half of the candidates will qualify. The first digit is already fixed by our 
selection of intervals. For the remaining k-1 digits we have ten possibilities, 
namely 0,1,2,3,4,5,6,7,8 and 9 of which only 2,3,5 and 7 are good. The 
probability that all k*l digits are good is therefore (4/10) kl . The probability q 
that a candidate qualifies as a member of the sequence is 



H 2 v io' 



3) 



The estimated number of members of the sequence in the interval (lO^lO** 1 ) 



is therefore given by m(k)=q-n(k) The estimated values are given in table 5. 
A comparison between the computer count and the theoretically estimated 
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values shows a very close fit as can be seen from diagram 5 where k>g 10 m is 
plotted against k. 



The prime-digital sub sequence 




k 



Diagram 5. Jogio m as a function of k. The upper curve corresponds to the 
computer count. 



loga logi 

For large values of k we can ignore the terms — 7 — and — 7 — in 

k k 

comparison with log 10 in (1). For large k we therefore have 



n(a,b,k) * 



(a -3)10* 
£logl0 



(!’) 



and ( 2 ) becomes 
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»(*)* 



4 10* 
A: log 10 



Combining this with (3) we get 



m(k ) * 



S-2 2 * 

JHoglO 



( 2 ’) 



(4) 



From which we see (apply for instance l'Hospital’s iuk) that m(k)— xx> as 
k-x». A fortiori the prime-digital sub sequence is infinite. 
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Chapter HI 



Non-Recursive Sequences 



UL1 Smadick Primitive Numbers 

Definition: For a given prime p and positive integer n the Smarandache 
primitive number Sp(n) is the smallest positive integer such that Sp(n>! is 
divisible by p“. 

This is important for the calculation of the Smarandache function 

S(n) which is defined as the wnallesi integer such that S(n)! is divisible by n. 
We note that S^n^S^®). 

It follows immediately from the definition that the sequence Sp(n) for a given 
prime p consists of multiples of p. Furthermore the definition implies that the 
sequence is non-descending, i.e. Sp(n)£ Sp(n+1). 

It is evident from the definition that 

Sp(n)=np for n£p (1) 

Upper and lower bounds for Sp(n) have been established by P51 GrftnSs [1 J: 

(p-l)n+l5 Sp(n)<pn (2) 

Let’s asssume Sp(n-l)<Sp(n) and p^l Sp(n) then it follows that Sp(n+k-l)= 
Sp(n+k-2)= ... = Sp(n+1)= Sp(n). This observation is expressed in the 
following way in Some Notions and Questions in Number Theory [2] 
Curious property: this is the sequence of multiples of p, each number being 
repeated as many times as its exponent (of power p) is.” 
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Algorithms for the calculation of Sp(n) are important as a stepping stone for 
the calculation of the Smarandache function. The lack of a closed formula 
makes this a very interesting topic. Three methods will be described, 
implemented and compared. 

Metbodl. 

Let’s denote Sp(n)=m. We want to calculate the smallest value of m for which 
1-2*3- ... -msO (mod p“). This is carried out in the following Ubasic program: 

10 input "p,n";P,N 
20 if N<=P then M 
30 Y-12=P A N 
40 while Y>0 
50 inc M 
60 Y=Y*M@Z 
70 wend 
80 print P,N,M 
90 end 

This program is effective for small values of n. As we move on to large 
exponents the increasingly large modulus slows down the program. A study of 
the compexity of this short and elegant program, which can equally well be 
used to calculate S(n), has been carried out by S and T. Tabirca [3], 

Method 2. 

When n<p we have Sp(n)=np. For n>p our potential solutions are multiples of 
p. The first multiple we need to examine is p 2 . For a given multiple m of p we 
determine first the largest power of p which is less than m, i.e. we determine 
k so that p k ^m<p lc+1 . We that continue by counting the number of factors p in 
m!. This is given by S = [m/p]+ [m/p 2 ] + [m/p 3 ] + ... + [m/p k ]. If S<n then wc 
proceed to the next multiple of p and so on until we arrive at S=0 in which 
case m=Sp(n). 

The following UBASIC program has been used to tabulate Sp(n) as well as 
the Smarandache function in the Smarandache Function Journal [4], 



=N*P:goto 80 



‘The factorial is reduced mod z each time 
'the loop is executed. 
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10 input "p,n";P,N 

20 if N<=P then M=N*P:goto 130 *Sp(n)=np in this case 
30 |=P-1 ‘i =p-1 so that m starts at p 2 

40 while S<N ‘in line 60 

50 inc I 

60 M=I*P ‘The starting point for m=p2 

70 S=0 

80 K=floor(kjg(M)/log(P)) 'Determine the largest k for which 
90 for J=1 to K ‘p^m 

100 S=S+floor(M/P A J) ‘Count and add the number times 

110 next 'pi occurs in m! 

120 wend 
130 print P,N,M 
140 end 

Method 3. 

This method is based on a theorem due to C. Dumitrescu and V.Seleacu {5], 
They prove the following: 

Theon*™ If for a given prime p the integer n is expressed in the number base 
bj=l, br=l+p, b3=l+p+p 2 , .... bk=l+p+ ... p kl , resulting in n=c,bi+ 0^+ ... 
+ cj\, where the first non-zero coefficient c^p and Cj<p for then 

Sp(n)=cip + C 2 P 2 + ... Ckp k . 

In the program below the number base is generated through the recursive 
formula bk=p-bk-i+l. 

10 dim C(100),B(100) 

20 input "p,n M ;P.N 
30 if N<=P then M=P*N:goto 150 
40 B(1)=1:l=l 
50 white N>=B(!) 

60 inc I 

70 B(I)=B(H)*P+1 'Calculate bz. t>3, etc 

80 wend 
90 K=l-1 

100 for l=K to 1 step -1 
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110 C(I}=N\B(I) 

120 N=res 
130 M=M+C(l)*PAi 
140 next 
150 print M 
160 end 



‘Calculate the coefficients 

‘Cl. C2, ... 

'Express Sp(n) using these 
'coefficients 



Remits. 

Which of these three programs is the most effective? Since Sp(n)=np when 
n^p all three programs will be equally efficient for n£p. The programs woe 
therefore tested by calculating Sp(n) for n=pt-l, p+2, ... p+50 for 
p=2,3,5,7,...229 (the first 50 primesj.The largest number manipulated in the 
programs is therefore 229 179 , which is a 559-digit number. The programs 
were fitted as routines in the same in- and output program. Short integers 
woe used whenever possible. To calculate these 2500 values of Sp(n) method 
I took 31 m 26 s, method II 22 s and method III II s on a Pentium 100 Mhz 
laptop. The rime t was measured in milliseconds. In figure 1 hi / is plotted 
against the number of primes for which the programs have been executed. It 
is seen that method in is the most effective. It runs about twice as fast as 
method II. Method I is effective fix' small primes but slows down considerably 
when the modulus p n increases. 



Int 




Figure 1. Comparison of execution times for methods I to III. 
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In table 1 consecutive squares have been emphasized by frames. The series 
light frames follow the rule Sp(p)= Sp(p+l)=p 2 which is obvious from 
definition. The series in bold frames is less obvious and this study will 
concluded by proving the following: 

Sp(4pt-4)= Sp(4jH-3)=4p 2 or S<p^’ ) )=S(p^ 3> )=(2p) 2 

Proof: Applying the theorem ofDumitrescu and Seleacu we have: 

4pM=4~(l+p), ie. Ci=0 and C2=4 from which 

S P (4p+p)=0-p+4-p 2 =(2p) 2 

and 

4p^3=p-l+3*(l+p), Le. Ci=p and Ci=3 from which 
S P (4p+3)=p-p+3-p 2 =(2p) 2 



III. 2 The Smarsuidache Function S(n) 

Definition: S(n) is the smallest integer such that S(n)! is divisible by n. 

The properties of this function have been subject to detailed studies in many 
papers. Some of these properties are listed here. When needed n is 
represented in the form n = pj 1 ' - pf 2 

1. S(n)£n 

2. p is prime iff S(p)=p 

3. If mn*0, then S(mn)sS(n)S(m) 
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4. If (m,n)=l, then S(mn)=max{S(m),S(n)} 

5. S(n)=niax{S Pi (a,), S Pi (a 2 ), ..\(a k )}, where Sy(n) are 

S marandacfae primitive numbers as defined in section II. 1 



6. If n is square free, Le. 04=02= ...=a*, then S(n)=max{pi, P2, ...pk} 

The purpose of this section, however, is to give an effective c omp u t e r 
algorithm 1 for the calc ulatio n of S(n) and to provide a table for die first 1000 
values of this function. 

10 'A U BASIC program for online calculation of S(n) 

20 dim F(10),E%(20) 

30 input "N";N 
40 gosub 70 
50 print “S(n)=“:S 
60 end 

70 N1=N:gosub 170 

N is factorized and the prime factors stored in F() with the 
corresponding exponents in E%() 

80 if F( 1 ) =N then S=N:goto 1 60 

tf F(l)=l then n is prime and consequently S(n)=n 
90 for J%=1 to K% 

100 if E%(J%)>1 then P2=F(J%):N2%=E%(J%):gosub 260:F{J%)=M 
1 10 next 

For prime factors of n with exponents larger than 1 F(J%)=P2 ts 
replaced by F(j%}=M which is calculated by one of the routines 
described eariier (here method 2) 

120 S=0 



1 The author was inspired to do this after several requests to help provide values for this function for 
Urge values of n. 
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130 for J%=1 to K% 

140 if F(J%)>S then S=F(J%) 

150 next 

The largest value ofF() is the value ofS(n) 

150 return 

1 70 ’Factorisation of N 1 , output f (i)Ae%(i) 

A simple routine to factorize N1 

180 K%=0rforl%=1 to 10:F(l%)=0:next 
190 forl%=1 to 20:E%(l%)=0:next 
200 while Nl>1 
210 P=prmdiv(Nl) 

220 if P=F(K%) then inc E%(K%) ebe inc K%:F(K%)=P:inc E%(K%) 
230 N1=N1\P 
240 wend 
250 return 

260 'Calculation of Sp(n). In: p2,n2%. Out: m. This routine is 
documented in section II. i 
270 I2=0:S2%=0 
280 while S2%<N2% 

290 inc 12 
300 M=I2*P2 
310 S2%=0 

320 K2=floor( log ( M) /log ( P2) } 

330 for J2=1 to K2 

340 S2%=S2%+floor(M/P2AJ2} 

350 next 
360 wend 
370 return 

This is the bare minimum program. When used in a special problem the 
number of instructions tend to increase considerably because of “safety 
valves” and data flow from and to other processes. 
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Table 2. The Smarandache function 



n 


0 


1 


_2 


_3 


_4 


5 


< 


7 


_8 


_8 






0 


2 


3 


4 


5 


3 


7 


4 


6 


1 


5 


11 


4 


13 


7 


5 


6 


17 


6 


19 


2 ~ 


5 


7 


11 


23 


4 


10 


13 


9 


7 


29 


3 


5 


31 


8 


11 


17 


7 


6 


37 


19 


13 


4 


5 


41 


7 


43 


11 


6 


23 


47 


6 


14 


5 


10 


17 


13 


53 


9 


11 


7 


19 


29 


59 


4 


5 


61 


31 


7 


8 


13 


11 


67 


17 


23 


7 


7 


71 


6 


73 


37 


10 


19 


11 


13 


79 


8 


6 


9 


41 


83 


7 


17 


43 


29 


11 


89 


9 


6 


13 


23 


31 


47 


19 


8 


97 


14 


11 


io_ 


10 


101 


17 


103 


13 


7 


53 


107 


9 


109 


11* 


11 


37 


7 


113 


19 


23 


29 


13 


59 


17 


12™" 


5 


22 


61 


41 


31 


15 


7 


127 


8 


43 


13 


13 


131 


11 


19 


67 


9 


17 


137 


23 


139 


14"™ 


7 


47 


71 


13 


6 


29 


73 


14 


37 


149 


15 


10 


151 


19 


17 


11 


31 


13 


157 


79 


53 


14 


8 


23 


9 


163 


41 


11 


83 


167 


7 


26 


17 


17 


19 


43 


173 


29 


10 


11 


59 


89 


179 


It 


€ 


181 


13 


61 


23 


37 


31 


17 


47 


9 


19 


19 


191 


8 


193 


97 


13 


14 


197 


11 


199 


20™" 


10 


67 


101 


29 


17 


41 


103 


23 


13 


19 


21 


7 


211 


53 


71 


107 


43 


9 


31 


109 


73 


22 


11 


17 


37 


223 


8 


10 


113 


227 


19 


229 


23_ 


23 


11 


29 


233 


13 


47 


59 


79 


17 


239 


24 


6 


241 


22 


12 


61 


14 


41 


19 


31 


83 


25 


15 


251 


7 


23 


127 


17 


10 


257 


43 


37 


24 


13 


29 


131 


263 


11 


53 


19 


89 


67 


269 


27"™ 


9 


271 


17 


13 


137 


11 


23 


277 


139 


31 


28 


7 


281 


47 


283 


71 


19 


13 


41 


8 


34 


29"“ 


29 


97 


73 


293 


14 


59 


37 


11 


149 


23 


30 


10 


43 


151 


101 


19 


61 


17 


307 


11 


103 


31 


31 


311 


13 


313 


157 


7 


79 


317 


53 


29 


32* 


8 


107 


23 


19 


9 


13 


163 


109 


41 


47 


33 


11 


331 


83 


37 


167 


67 


7 


337 


26 


113 


34 


17 


31 


19 


21 


43 


23 


173 


347 


29 


34 9 


35“ 


10 


13 


11 


353 


59 


71 


89 


17 


179 


359 


34 


6 


38 


181 


22 


13 


73 


61 


367 


23 


41 


37™" 


37 


53 


31 


373 


17 


15 


47 


29 


9 


379 


38 


19 


127 


191 


383 


8 


11 


193 


43 


97 


389 


39 


13 


23 


14 


131 


197 


79 


11 


397 


199 


19 


40 


10 


401 


67 


31 


101 


9 


29 


37 


17 


409 


41 


41 


137 


103 


59 


23 


83 


13 


139 


19 


419 


42 


7 


421 


211 


47 


53 


17 


71 


61 


107 


13 


43 


43 


431 


9 


433 


31 


29 


109 


23 


73 


439 


44 


11 


14 


17 


443 


37 


89 


223 


149 


8 


449 


45™ 


10 


41 


113 


151 


227 


13 


19 


457 


229 


17 


44 


23 


461 


11 


463 


29 


31 


233 


4 67 


13 


67 


47 


47 


157 


59 


43 


79 


19 


17 


53 


239 


479 


48 


8 


37 


241 


23 


22 


97 


12 


487 


61 


163 


49 


14 


491 


41 


29 


19 


11 


31 


71 


83 


499 
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SO ! 15 167 251 503 7 101 23 26 127 509 

Table 2 ctd. The Smarandache function 



5 « 7 8 » 



51 


17 


73 


12 


19 


257 


103 


43 


47 


37 


173 


52 


13 


521 


29 


523 


131 


10 


263 


31 


11 


46 


53 


53 


59 


19 


41 


89 


107 


67 


179 


269 


14 


54 


9 


541 


271 


181 


17 


109 


13 


547 


137 


61 


55 


11 


29 


23 


79 


277 


37 


139 


557 


31 


43 


56 


7 


17 


281 


563 


47 


113 


283 


9 


71 


569 


57 


19 


571 


13 


191 


41 


23 


8 


577 


34 


193 


58 


29 


83 


97 


53 


73 


13 


293 


587 


14 


31 


58 


59 


197 


37 


593 


11 


17 


149 


199 


23 


599 


•O 


10 


601 


43 


67 


151 


22 


101 


607 


19 


29 


Cl 


61 


47 


17 


613 


307 


41 


11 


617 


103 


619 


82 


31 


23 


311 


89 


13 


20 


313 


19 


157 


37 


83 


7 


631 


79 


211 


317 


127 


53 


14 


29 


71 


64 


8 


641 


107 


643 


23 


43 


19 


647 


9 


59 


65 


13 


31 


163 


653 


109 


131 


41 


73 


47 


659 


66 


11 


661 


331 


17 


83 


19 


37 


29 


167 


223 


67 


67 


61 


8 


673 


337 


10 


26 


€77 


113 


97 


68 


17 


227 


31 


683 


19 


137 


21 


229 


43 


53 


69 


23 


691 


173 


11 


347 


139 


29 


41 


349 


233 


70 


10 


701 


13 


37 


11 


47 


353 


101 


59 


709 


71 


71 


79 


89 


31 


17 


13 


179 


239 


359 


719 


72 


6 


103 


38 


241 


181 


29 


22 


727 


13 


15 


73 


73 


43 


61 


733 


367 


14 


23 


67 


41 


739 


74 


37 


19 


53 


743 


31 


149 


373 


83 


17 


107 


75 


15 


751 


47 


251 


29 


151 


9 


757 


379 


23 


76 


19 


761 


127 


109 


191 


17 


383 


59 


10 


769 


77 


11 


257 


193 


773 


43 


31 


97 


37 


389 


41 


78 


13 


71 


23 


29 


14 


157 


131 


787 


197 


263 


79 


79 


113 


11 


61 


397 


53 


199 


797 


19 


47 


80 


10 


89 


401 


73 


67 


23 


31 


269 


101 


809 


81~ 


9 


811 


29 


271 


37 


163 


17 


43 


409 


13 


82~ 


41 


821 


137 


823 


103 


11 


59 


827 


23 


829 


83 


83 


277 


13 


17 


139 


167 


19 


31 


419 


839 


84 


7 


58 


421 


281 


211 


26 


47 


22 


53 


283 


85 


17 


37 


71 


853 


61 


19 


107 


857 


13 


859 


86 


43 


41 


4 31 


863 


9 


173 


433 


34 


31 


79 


87 


29 


67 


109 


97 


23 


15 


73 


877 


439 


293 


88 


11 


881 


14 


883 


17 


59 


443 


887 


37 


127 


89 


89 


11 


223 


47 


149 


179 


8 


23 


449 


31 


90 


10 


53 


41 


43 


113 


181 


151 


907 


227 


101 


91 


13 


911 


19 


83 


457 


61 


229 


131 


17 


919 


92 


23 


307 


461 


71 


11 


37 


463 


103 


29 


929 


93 


31 


19 


233 


311 


467 


17 


13 


937 


67 


313 


94 


47 


941 


157 


41 


59 


9 


43 


947 


79 


73 


95 


19 


317 


17 


953 


53 


191 


239 


29 


479 


137 


96 


8 


62 


37 


107 


241 


193 


23 


967 


22 


19 


97 


97 


971 


12 


139 


487 


13 


61 


977 


163 


89 


98 


14 


109 


491 


983 


41 


197 


29 


47 


19 


43 



57 






H : 11 991 31 331 71 199 83 997 499 37 

1IL3 Smirudack m-Power Residues 

Definition: The m-power residue of n is the largest m-power free number 
which divides n. Notation H<n). 



Let’s express n in the form n = pf* -pf 2 - -Pi‘ then M r (n) = p“' p® 1 -...p“ k , 
where a 1 =min{a,,m-l} for i=l, 2, ... k. 

It follows directly from the definition that KL(n) is a multiplicative arithmetic 
function, which is expressed in (I) below. 

If (n,m)=l then M^nm^M^nl-M^m) 

We will only make calc ulatio ns for two special cases: Cubical residues and 
Square residues. 



IEOa Cubical Residues. 

In the prime factor representation of n a£3 will be replaced by a,— 2. A 
simple Ubasic routine results in table 3. 



Table 3. The first 100 terms of the cubical residues sequence. 



n 


0 


1 


2 


3 


4 


5 


6 


7 


8 


9 






1 


2 


3 


4 


5 


6 


7 


4 


9 


1 


10 


11 


12 


13 


14 


15 


4 


17 


18 


19 


2 


20 


21 


22 


23 


12 


25 


26 


9 


28 


29 


3 


30 


31 


4 


33 


34 


35 


36 


37 


38 


39 


4 


20 


41 


42 


43 


44 


45 


46 


47 


12 


49 


5 


50 


51 


52 


53 


18 


55 


28 


57 


58 


59 


6 


60 


61 


62 


63 


4 


65 


66 


67 


68 


69 


7 


70 


71 


36 


73 


74 


75 


76 


77 


78 


79 


8 


20 


9 


82 


83 


84 


85 


86 


87 


44 


89 


9 


90 


91 


92 


93 


94 


95 


12 


97 


98 


99 



ss 







m~3b Square Residues 



For n in the form n = p“‘ p? 2 -...p^‘ foe square residue SJn) is given by 

S^n)=p,.p2- ... |\. 

Apart from foe multiplictivc property we add two more pr op er ties of this 
sequence. 

If aX) and p any prime number then Sr<p a )=p. 

SXnm^SX^-SXm) /(Sr(n), SXm)). 

The first 100 terms af this sequence are given in table 4. 



Table 4 . The first 100 terms of the square residues sequence 



n 


0 


1 


2 


3 




5 


6 


7 


8 


9 






1 


2 


3 


2 


5 


6 


7 


2 


3 


1 


10 


11 


6 


13 


14 


15 


2 


17 


6 


19 


2 


10 


21 


22 


23 


6 


5 


26 


3 


14 


29 


3 


30 


31 


2 


33 


34 


35 


6 


37 


38 


39 


4 


10 


41 


42 


43 


22 


15 


46 


47 


6 


7 


5 


10 


51 


26 


53 


6 


55 


14 


57 


58 


59 


6 


30 


61 


62 


21 


2 


65 


66 


67 


34 


69 


7 


70 


71 


6 


73 


74 


15 


38 


77 


78 


79 


8 


10 


3 


82 


83 


42 


85 


86 


87 


22 


89 


9 


30 


91 


46 


93 


94 


95 


6 


97 


14 


33 



In chapter VI we will return to the Smarandache function and the square 
residues. 
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59 







Chapter IV 



Periodic Sequences 



IV.l Introduction 

In Mathematical Spectrum, vd 29 No 1 [1], is an article on Smarandache’s 
periodic sequences which terminates with the statement. 

“There will always be a periodic sequence whenever we have a 
function f:S->S, where S is a finite set of positive integers and we 
repeat the function f. ” 

However, consider the following trivial function f(Xk):S-»S, where S is an 
ascending set of integers fo, 32, ... a,, ... a„}: 

f Xu if 

f(Xk)=H x* if Xi=a, 
l Xfc+i ifxk<a. 

As we can see the iteration of the function f in this case converges to an 
invariant a^ which we may of course consider as a sequence (or loop) of only 
one member. In this study, however, a periodic sequence will be referred to 
as if it has at least two members. If it has only one member it will be 
referred to as invariant 

There is one more saag to overcome. In the Smarandache sequences OS is 
considered as a two-digit integer. The consequence of this is that 00056 is 
considered as as a five digit integer while 056 is considered as a three-digit 
integer. We will abolish this ambiguity, 05 is a one-digit integer and 00200 is 
a three-digit integer. 
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With these two remarks in mind let’s look at these sequences. There are in all 
four different ones reported in the above mentioned article in Mathematical 
Spectrum. The study of the first one will be carried out in much detail in view 
of the above remarks. 



IV.2 The Two-Digit Smarandache Periodic Sequence 

It has been assumed that the definition given below leads to a repetition 
according to Dirichlet’s box principle (or the statement made above). 
However, as we will see, this definition leads to a collapse of the sequence. 

Preliminary definition. Let N k be an integer of at most two digits and let N k ’ 
be its digital reverse. We define the element of the sequence through 



Nfcn= I N k -N k ’ | 



where the sequence is initiated by an arbitrary two digit integer N] . 

Let’s write Ni in the form Ni=10a+b where a and b are digits. We then have 
N 2 = I lOa+b-lOb-a | =9- 1 a-b | 



The | a-b | can only assume 10 different values 0,1,2, ... ,9. This means that 
N 3 is generated from only 10 different values of N 2 . Let’s first find out which 
two digit integers result in | a-b I = 0,1,2, .. and 9 respectively. 



I a-b | Corresponding two digit integers 



0 


11 


22 


33 


44 


55 


66 


77 


88 


99 


1 


10 


12 


21 


23 


32 


34 


43 


45 


54 


2 


13 


20 


24 


31 


35 


42 


46 


53 


57 


3 


14 


25 


30 


36 


41 


47 


52 


58 


63 


4 


15 


26 


37 


40 


48 


51 


59 


62 


73 


5 


16 


27 


38 


49 


50 


61 


72 


83 


94 


6 


17 


28 


39 


60 


71 


82 


93 






7 


19 


29 


70 


81 


92 










8 


19 


80 


91 














9 


90 



















56 


65 


67 


76 


78 


87 89 98 


64 


68 


75 


79 


86 


97 


69 


74 


85 


96 






84 


95 
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It is now easy to follow the iteration of the sequence which invariably 
terminates inO, table 1. 

1 - Iteration of secjjence accofdr>g to the preiminary detrition 



■BUI 


Nj 


Ns 


N 4 


Ns 


Ns 


Ns 


0 


0 












1 


9 


0 










2 


18 


63 


27 


45 


9 


0 


3 


27 


45 


9 


0 






4 


36 


27 


45 


9 


0 




5 


45 


9 


0 








6 


54 


9 


0 








7 


63 


27 


45 


9 


0 




8 


72 


45 


9 


0 






9 


81 


63 


27 


45 


9 


0 



The termination of the sequence is preceded by the one digit element 9 whose 
reverse is 9. The following definition is therefore proposed. 

Definition of Smarandachc’s two-digit periodic sequence. Let N k be an 
integer of at most two digits. N k ’ is defined through 



Ithe reverse of Nk if N k is a two digit integer 

Nk’ = *{ 

lNk‘10 if N k is a one digit integer 
We define the element Nk+i of the sequence through 
Nk+i = | NfNk’ I 

where the sequence is initiated by an arbitrary two digit integer N, with 
unequal digits. 

Modifying table 1 according to the above definition results in table 2. 
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Conclusion: The iteration always produces a loop of length 5 which starts on 
the second or the third term of the sequence. The period is 9, 81, 63, 27, 45 or 
a cyclic permutation thereof. 



Table 2. Iteration of the Smarandache two digit sequence 



a-b| 


Nz 


Ns 


N4 


Ns 


Ns 


Ns 


N 7 


1 


9 


81 


63 


27 


45 


9 




2 


18 


63 


27 


45 


9 


81 


63 


3 


27 


45 


9 


81 


63 


27 




4 


36 


27 


45 


9 


81 


63 


27 


5 


45 


9 


81 


63 


27 


45 




6 


54 


9 


81 


63 


27 


45 


9 


7 


63 


27 


45 


9 


81 


63 




8 


72 


45 


9 


81 


63 


27 


45 


9 


81 


63 


27 


45 


9 


81 





IV.3 The Smarandache n-digit periodic sequence. 

Let’s extend the definition of the two-digit periodic sequence in the following 
way. 

Definition of the Smarandache n-digit periodic sequence. 

Let Nk be an integer of at most n digits and let R* be its reverse. Nk’ is 
defined through 

N' k =R k - io'-'T 10810 ^) 



We define the element of the sequence through 



N w =|N k -N k ’| 



where the sequence is initiated by an arbitrary n-digit integer Nj in the 
domain 10 n ^N 1 <10 It ^ 1 . It is obvious from the definition that O^N^IO^ 1 , 
which is the range of the iterating function. 
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Let’s consider the cases n=3, n=4, n=5 and n=6. 

n*3. 

Domain 100^Ni£999. . The iteration will lead to an invariant or a loop 
(periodic sequence) 1 . There are 90 symmetric integers in the domain, 101, 
111, 121, ...202, 212, ..., for which Nj=0 (invariant). All other initial 
integers iterate into various entry points of the same periodic sequence. The 
number of numb ers in the domain resulting in each entry of the loop is 
denoted s in table 3. 

Table 3. Smarandache 3-dgit periotic sequence 

s 239 H 200 240 ] 20 ~ 

Loop 99 891 693 297 495 



It is easy to explain die relation between this loop and the loop found for n=2. 
Consider N=ao+10ai+100a2. From this we have [N-N’|=99|a2-ao|=I l-9|a2-ao| 
which is 11 times the corresponding expression for n=2 and as we can see this 
produces a 9 as middle (or first) digit in the sequence for n=3. 



Domain 1000^NiS9999. The largest number of iterations carried out in order 
to reach the first member of the loop is 18 and it happened for N t =1019. The 
iteration process ended up in the invariant 0 for 182 values of Nj, 90 of these 
are simply the symmetric integers in the domain like Ni=4334, 1881, 7777, 
etc., the other 92 are due to symmetric integers obtained after a couple 
iterations. Iterations of the other 8818 integers in the domain result in one 
the following 4 loops or a cyclic permutation of one of these. The number 
numbers in the domain resulting in each entry of the loops is denoted s 
table 4 



1 TIb* is dabonled a 4n»il m Surfing an At Ocean ofNumben by the author, V«il Univ. Pre*» 
1997. 
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Table 4. Smarandache 4-digit periodic sequences 



s 

Loop 


378 

2178 


259 

6534 








s 


324 


18 


288 


2430 


310 


Loop 


90 


810 


630 


270 


450 


s 


446 


2 


449 


333 


208 


Loop 


909 


8181 


6363 


2727 


4545 


s 


329 


11 


290 


2432 


311 


Loop 


999 


8991 


6993 


2997 


4995 



n=5. 



Domain 10000^Ni^99999. There are 900 symmetric integers in the domain. 
920 integers in the domain iterate into the invariant 0 due to symmetries. 

Table 5. Smarandache 5-digit periodic sequences 



s 

Loop 


3780 

21978 


2590 

65934 








s 


3240 


180 


2880 


24300 


3100 


Loop 


990 


8910 


6930 


2970 


4950 


s 


4469 


11 


4490 


3330 


2080 


Loop 


9009 


81081 


63063 


27027 


45045 


s 


3299 


101 


2900 


24320 


3110 


Loop 


9999 


89991 


69993 


29997 


49995 



n=6. 

Domain 100000^^999999. There are 900 symmetric integers in the 
domain. 12767 integers in the domain iterate into the invariant 0 due to 
symmetries. The longest sequence of iterations before arriving at the first loop 
member is 53 for N= 100720. The last column in table 6 shows the number of 
integers iterating into each loop. 
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Table 6. Smarandache 6-cSgit periodic sequences 



s 

L 1 


13667 

0 


13667J 


s 


13542 


12551 


















L 2 


13586 


65340 
















■HI 


s 


12685 


12685 
















26271 


L 3 


219978 


659934 


















n 


19107 


2711 


7127 


123320 


12446 










164711 


H 


900 


8100 


6300 


2700 


4500 












n 


25057 


18 


12259 


20993 


4449 










62776 


Ifl 


9090 


81810 




27270 


45450 










1H 


o 


47931 


14799 


42603 


222941 


29995 










3582691 


m 


9990 


89910 


69930 


29970 


49950 













a 


25375 


11 


12375 


21266 


4409 












119 


90009 


810061 


630063 


270027 


450045 










■HI 




1488 


2 


1005 


1033 


237 










3765 


18 


90909 


818181 


636363 


272727 


454545 










■■ ■ 




1809 


11 


1350 


1570 


510 














99099 


891891 


693693 


297297 


495495 










■■ 




19139 


2648 


7292 


123673 


12472 












L 10 


99999 


899991 


699993 


299997 


499995 










■ ■ 


s 


152 


4 


1254 


972 


492 


in 


826 


485 


429 


4725 


LI) 


10989 


978021 


857142 


615384 


131868 


736263 


373626 


252747 


494505 


L—J 


s 


623 


64 


156 


796 


377 


36 


525 


140 


194 


■■■ 


L 12 


43659 


912681 


726462 


461835 


76329 


847341 


703593 


306286 


374517 


■ 


ctd 


596 


117 


156 


793 


327 


65 


530 


139 


179 


58131 


c\ d 


340956 


318087 


462726 


164538 


670923 


341847 


406296 


286306 


517374 


HU 



Table 7. The subtraction periodic sequence. 10<Ni599 



*1 


12 


20 


1 


9 


89 


97 


78 


86 


67 


75 


56 


64 


45 


53 


34 


42 


23 


31 


*2 


13 


30 


2 


19 


90 


8 


79 


96 


68 


85 


57 


74 


46 


63 


35 


52 


24 


41 


*3 


14 


40 


3 


29 


91 


18 


80 


7 


69 


95 


58 


84 


47 


73 


35 


62 


25 


51 


*4 


15 


50 


4 


39 


92 


28 


81 


17 


70 


6 


59 


94 


48 


83 


37 


72 


26 


61 


*5 


16 


60 


5 


49 


93 


38 


82 


27 


71 





















IV. 4 The Sairudxhc Subtraction Periodic Sequence 

Definition: Let N k be a positive integer of at most n digits and let Rt be its 
digital reverse. Nk’ is defined through 
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N k = R k ■\0’ t -'~ [ ' oi>tNl] 

We define the element N^, erf the sequence through 

Nw=|tV-cl 

where c is a positive integer. The sequence is initialed by an aibitraiy positive 
n -djgit integer Nj. It is obvious from the definition that 0^N k <10^ , which is 
the range of the iterating function. 

c-l,«=2,l©i£Ni£99 

When N, is of the form Ilk or llk-1 then the iteration process results in 0, 
see figure la. 

Every other member of the interval 1<KN,S99 is a entry point into one of five 
Hifferrnt cyclic periodic sequences. Four of these are of length 18 and one of 
length 9 as shown in table 7 and illustrated in figures lb and lc, where 
important features of the iteration chains are shown. 
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I5c^9, n=2, 100£N,iJ999 



A computer analysis revealed a number of interesting facts concerning the 
application of the iterative function. 

There are no periodic sequences for c=l, c=2 and c=5. All iterations result in 
the invariant 0 after, sometimes, a large number of iterations. 

Table 8. Loop statistics. L=4ength of loop, f=first term of loop 
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For the other values of c there are always some values of N, which do not 
produce periodic sequences but terminate on 0 instead. Those values of N! 
which produce periodic sequences will either have N] as the first term of the 
sequence or one of the values f determined by lsfsc-l as first term. There are 
only eight different possible value for the length of the loops, name ly 12 , 22, 
33, 50, 100, 167, 189, 200. Table 8 shows how many of the 900 initi atin g 
integers in the interval 100^N,^999 result in each type of loop or invariant 0 
for each value of c. 



A few examples: 

For c=2 and N,=202 the sequence ends in the invariant 0 after only 2 
iterations: 

202 200 0 

For c=9 and N]=208 a loop is closed after only 1 1 iterations: 

208 793 388 874 469 955 550 46 631 127 712 208 
For c=7 and N]=109 we have an example of the longest loop obtained. It has 
200 elements and the loop is closed after 286 iterations: 



109 894491 187 774 470 67 753 350 46 633 329 916 612 209 895 591 188 874471 
167 754 450 47 733 330 26 613 309 896 691 189 974 472 267 755 550 48 833 331 
126 614 409 897 791 190 84 473 367 756 650 49 933 332 226 615 509 898 891 191 
184 474 467 757 750 50 43 333 326 616 609 899 991 192 284 475 567 758 850 51 
143 334 426 617 709 900 2 1 93 384 476 667 759 950 52 243 335 526 618 809 901 
1 02 1 94 4§4 4 77 767 760 60 53 343 336 626 619 909 902 202 195 584 478 867 761 
160 54 443337 726 620 19 90 3 302 196 684 479 967 762 260 55 543 338 826 621 
1 19 9Q4 4Q2 197 784 480 77 763 360 56 643 339 926 622 219 905 SO? IPS wa 4a] 
1 77 764 46Q 57 743 340 36 623 319 9 06 602 1 99 984 482 277 765 560 58 843 34] 
136 624 41 9 907 702 200 5 493 387 776 670 69 953 352 746 635 529 91 3 SI ? 91 1 
1 05 494 487 777 770 70 63 35 3 346 636 629 919 912 212 205 495 587 778 870 71 
Ii3 354 446 637 729 920 2? 21 3 305 496 687 779 970 72 263 355 546 638 829 921 
1 22 21 4 405 497 787 780 80 73 363 356 646 639 929 92? ??? 91 5, SOS 498 887 781 
1 80 74 463 35 7 746 6 4Q 3 9 923 322 21 6 605 499 987 782 280 75 563 358 846 64 1 
139 924 422 217 705 500 2 



IV.5 The Smarandache Multiplication Periodic Sequence 

Definition: Let Ol be a fixed integer and No and arbitrary positive integer. 
Ntn is derived from N k by multiplying each digit x of N k by c retaining only 
the last digit of the product cx to become the corresponding digit of N k+1 . 
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In this each digit position goes through a separate development without 
interference with the surrounding digits. Let’s as an example consider the 
third digit of a 6-digit integer for c=3. The iteration of the third digit follows 
the schema - 

xxlyyy the third digit has been arbitrarily chosen to be 7. 

xxlyyy 

xx3yyy 

™9yyy 

xxlyyy which closes the loop for the third digit 

Let’s now consider all the digits of a six-digit integer 237456: 

237456 

691258 

873654 

419852 

237456 which closes the loop. 

The digits 5 and 0 are invariant under this iteration. All other digits have a 
period of 4 for c=3. 

Conclusion: Integers whose digits are all equal to 5 are invariant under the 
given operation. All other integers iterate into a loop of length 4. 

We have seen that the iteration process for each digit for a given value of c 
completely determines the iteration process for any n-digit integer. It is 
therefore of interest to see these angle digit iteration sequences: 

With the help of table 9 it is now easy to characterize the iteration process for 
each value of c. 

Integers composed of the digit 5 result in an invariant after one iteration. 
Apart form this we have for 

c=2. Four term loops starting on the first or second term. 
c=3. Four term loops starting with the first term. 

c=4. Two tram loops starting on the first or second term (could be called a 
switch or pendulum). 
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Table 9. One-digit multiplication sequences 



12486211 397 1 



2 4 8 6 2 



2 6 8 4 2 



362486 397 13 



4 8 6 2 4 

5 0 0 

6 2 4 8 6 



4 2 6 8 4 

5 5 

6 8 4 2 6 



748624 71 397 



8 6 2 4 8 



8 4 2 6 8 



98624819 7 139 



14 6 4 

2 8 2 

3 2 8 2 

4 6 4 

5 0 0 

6 4 6 

7 8 2 8 

8 2 8 

9 6 4 



1 5 5 

2 0 0 

3 5 5 

4 0 0 

5 5 

6 0 0 

7 5 5 

8 0 0 



1 6 6 
2 2 

3 8 8 

4 4 

5 0 0 

6 6 

7 2 2 

8 8 

9 4 4 



17 9 3 1 

2 4 8 6 2 

3 17 9 3 

4 8 6 2 4 

5 5 

6 2 4 8 6 

7 9 3 17 

8 6 2 4 8 

9 3 17 9 







c* 


=8 






1 


8 


4 


2 


6 


8 


2 


6 


8 


4 


2 




3 


4 


2 


6 


8 


4 


4 


2 


6 


8 


4 




5 


0 


0 








6 


8 


4 


2 


6 




7 


6 


8 


4 


2 


6 


8 


4 


2 


6 


8 




9 


2 


6 


8 


4 


2 



2 8 2 

3 7 3 

4 6 4 

5 5 

6 4 6 

7 3 7 

8 2 8 
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c=S. Invariant after (me iteration. 



c=€. Invariant after one iteration, 

c=7. Four term loop starting with the first term. 

c=8. Four term loop starting with the second term 

c=9. Two term loops starting with the first tram (pendulum). 



IV.6 The Smarandache Mixed Composition Periodic Sequence 

Definition. Let No be a two-digit integer arlO+ao. If ai+ao<10 then bi= ai+ao 
otherwise bi= ai+ao+1. bo=(ai-aoj . We define Ni=bi lCHbo. is derived 
from Nk in the same way. 2 



There are no invariants in this case. 36, 90, 93 and 99 produce two-element 
loops. The longest loops have 18 elements. A complete list of these periodic 
sequences is presented below. 



10 


n 


20 


22 


40 


44 


80 


88 


70 


77 


50 


55 


10 
















11 


20 


22 


40 


44 


80 


88 


70 


77 


50 


55 


10 


11 
















12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 




13 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


u 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 




15 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 




17 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


34 


34 


71 


86 


18 


97 


72 


95 


54 


91 


18 




























1? 


18 


97 


72 


95 


54 


91 


18 


























20 


22 


40 


44 


80 


88 


70 


77 


50 


55 


10 


11 


20 
















21 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


22 


40 


44 


80 


88 


70 


77 


50 


55 


10 


11 


20 


22 
















23 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


24 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 



2 Formulation conveyed to the author “Let N be* two-digit number. Add the digits, and add them 
again if the sum is greater than 10. Abo take the Absolute value of their difference. These are the first 
and second digits ofNi/* 
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25 73 14 53 82 16 75 

26 84 34 71 86 52 73 

27 95 54 91 18 97 72 

28 16 75 32 51 64 12 

29 27 95 54 91 18 97 

30 33 60 66 30 

31 42 62 84 34 71 86 

32 51 64 12 31 42 62 

33 60 66 30 33 

34 71 86 52 73 14 53 

35 82 16 75 32 51 64 

36 93 36 

37 14 53 82 16 75 32 

38 25 73 14 53 82 16 

39 36 93 36 

40 44 80 88 70 77 50 

41 53 82 16 75 32 51 

42 62 84 34 71 86 52 

43 71 86 52 73 14 53 

44 80 88 70 77 50 55 

45 91 18 97 72 95 54 

46 12 31 42 62 84 34 

47 23 51 64 12 31 42 

48 34 71 86 52 73 14 

49 45 91 18 97 72 95 

50 55 10 11 20 22 40 

51 64 12 31 42 62 84 

52 73 14 53 82 16 75 

53 82 16 75 32 51 64 

54 91 18 97 72 95 54 

55 10 11 20 22 40 44 

56 21 31 42 62 84 34 

57 32 51 64 12 31 42 

58 43 71 86 52 73 14 

59 54 91 18 97 72 95 

60 66 30 33 60 

61 75 32 51 64 12 31 

62 84 34 71 86 52 73 

63 93 36 93 

64 12 31 42 62 84 34 

65 21 31 42 62 84 34 

66 30 33 60 66 

67 41 53 82 16 75 32 

68 52 73 14 53 82 16 

69 63 93 36 93 



32 51 64 12 31 42 62 
14 53 82 16 75 32 51 
95 

31 42 62 84 34 71 86 

72 95 

52 73 14 53 82 16 75 
84 34 71 86 52 73 14 

82 16 75 32 51 64 12 
12 31 42 62 84 34 71 

51 64 12 31 42 62 84 
75 32 51 64 12 31 42 

55 10 11 20 22 40 
64 12 31 42 62 84 34 

73 14 53 82 16 75 32 
82 16 75 32 51 64 12 
10 11 20 22 40 44 

91 

71 86 52 73 14 53 82 
62 84 34 71 86 52 73 

53 82 16 75 32 51 64 

54 91 

44 80 88 70 77 50 
34 71 86 52 73 14 53 

32 51 64 12 31 42 62 
12 31 42 62 84 34 71 

80 88 70 77 50 55 
71 86 52 73 14 53 82 
62 84 34 71 86 52 73 

53 82 16 75 32 51 64 

54 

42 62 84 34 71 86 52 
14 53 82 16 75 32 51 

71 86 52 73 14 53 82 
71 86 52 73 14 53 82 

51 64 12 31 42 62 84 
75 32 51 64 12 31 42 



84 34 71 86 52 73 
64 12 31 42 62 84 

52 73 14 53 82 16 



32 51 64 12 31 

53 82 16 75 32 

31 42 62 84 34 

86 52 73 14 53 82 

34 71 86 52 73 14 

62 84 34 71 86 52 73 



71 86 52 73 14 53 
51 64 12 31 42 
31 42 62 84 34 71 



16 75 32 51 64 12 
14 53 82 16 75 32 51 
12 31 42 62 84 34 



82 16 75 32 51 
84 34 71 86 52 
86 52 73 14 53 



16 75 32 51 64 12 31 
14 53 82 16 75 32 
12 31 42 62 84 34 71 



73 14 53 82 16 75 

64 12 31 42 62 

16 75 32 51 64 

16 75 32 51 64 12 31 

34 71 86 52 73 14 53 

62 84 34 71 86 52 



73 




70 


77 


50 


55 


10 


11 


20 


22 


40 


44 


80 


88 


70 


















71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 






72 


95 


54 


91 


18 


97 


72 






























73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 






74 


23 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 






76 


41 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


77 


50 


55 


10 


11 


20 


22 


40 


44 


80 


88 


70 


77 


















78 


61 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


79 


72 


95 


54 


91 


18 


97 


72 




























80 


88 


70 


77 


50 


55 


10 


11 


20 


22 


40 


44 


80 


















81 


97 


72 


95 


54 


91 


18 


97 




























82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 






83 


25 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 






85 


43 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 






87 


61 


75 


32 


51 


64 


12 


31 


42 


62 


84 


34 


71 


86 


52 


73 


14 


53 


82 


16 


75 


88 


70 


77 


50 


55 


10 


11 


20 


22 


40 


44 


80 


88 


















89 


81 


97 


72 


95 


54 


91 


18 


97 


























90 


99 


90 






































91 


18 


97 


72 


95 


54 


91 






























92 


27 


95 


54 


91 


18 


97 


72 


95 


























93 


36 


93 






































94 


45 


91 


18 


97 


72 


95 


54 


91 


























95 


54 


91 


18 


97 


72 


95 






























96 


63 


93 


36 


93 


































97 


72 


95 


54 


91 


18 


97 






























98 


81 


97 


72 


95 


54 


91 


18 


97 


























99 


90 


99 
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Chapter V 

Smarandache Concatenated Sequences 



V.l Introduction 

Smarandache formulated a series of very artificially conceived sequences 
through concatenation. The sequences studied below are special cases of the 
Smarandache Concatenated S -sequence. 



Definition: Let G-{gi, g 2 , .... .... } be an ordered set of positive integers 

with a given property G. The corresponding concatenated S.G sequence is 
defined through 

S.G = {a,:a, = g x ,a k = a k _ x • io 1+log '°** + g k ,k> 1} 

In table 1 the first 20 terms are listed for three cases, which we will deal with 
in some detail below. 



V.2 The Smarandache Odd Sequence 

The Smarandache Odd Sequence is generated by choosing 

G = {1,3,5,7,9,11 }. Smarandache asks how many terms in this sequence 

are primes and as is often the case we have no answer. But for this and the 
other concatenated sequences we can take a look at a fairly large number of 
terms and see how frequently we find primes or potential primes. As in the 
case of prime-product sequence we will resort to Fermat’s little theorem to 
find all primes/pseudo-primes among the first 200 terms. If they are not too 
big we can then proceed to test if they are primes. For the Smarandache Odd 
Sequence there are only five cases which all were confirmed to be primes 
using the elliptic curve prime factorization program. In table 2 # is the term 
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number, L is the numb er of digits of N and N is a prime number member of 
the Smarandache Odd Sequence : 

Table la. The first 20 terms of the Smarandache Odd Sequence 



1 

13 
135 
1357 
13579 
1357911 
135791113 
13579111315 
1357911131517 
135791113151719 
13579111315171921 
1357911131517192123 
135791113151719212325 
13579111315171921232527 
1357911131517192123252729 
1357911 13151719212325272931 
13579111315171921232527293133 
1 3579 1 1 1 3 1 5 1 7 1 92 1 232527293 1 3335 
1 3579 1 1 1 3 1 5 1 7 1 92 1 232527293 1 333537 
1 3579 111315171 9212325272931 33353739 
1 3579 111315171921 232527293 1 333537394 1 



Table 2. Prime numbers in the Smarandache Odd SeQjence 



# 


L 


N 


2 


2 


13 


10 


15 


135791113151719 


16 


27 


135791 113151 71 921 2325272931 


34 


63 


135791 1 131517192123252729313335373941434547495153555 
7596)636567 


49 


93 


135791 1 131517192123252729313335373941434547495153555 
75961 6365676971 7375777981 8385878991 939597 
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Term #201 is a 548 digit number. 

V.3 The Snurandache Eva S eq uenc e 

The Snurandache Eves Sequence is generated by choosing G={2,4,6,8,10, 

}. The question here is : How many terms are nth powers of a positive 

integer? 

A tom which is a nth power must be of the form 2°-a where a is an odd n th 
power. The first step is therefore to find the highest power of 2 which divides 
a given member of the sequence, i.e. to determine n and at the nme time we 
will find a. We then have to test if a is a nth power. The Ubasic program 
below has been implemented for the first 200 toms of the sequence. No nth 
powers were found. 

Table lb. The first 20 terms of the Smarandache Even Sequence 



2 
24 
246 
2468 
246810 
24681012 
2468101214 
246810121416 
24681012141618 
2468101214161820 
246810121416182022 
24681012141618202224 
2468 101214161 820222426 
246810121416182022242628 
24681012141618202224262830 
2468101214161820222426283032 
246810121416182022242628303234 
24681012141618202224262830323436 
2468101214161820222426283032343638 



77 







2468 101214161 82022242628303234363840 
24681012141618202224262830323436384042 



Ubasic program: forty the essential part of the program is fated) 

60 N=2 

70 for U%*4 to 400 step 2 
80 D%=intOog(U%)/log(10))+l 
90 N=N*10AD%+U% 

100 A=N:E%=0 
110 repeat 

120 Al=A*=A\25nc E% 

IX until resoO 
132 dec E%l*=A 1 
140 B=round(AAfi/E%)) 

IX if &ae%= a then print E%.N 
160 next 
170 end 

So there is not even a perfect square among the first 200 terms of the 
Smarandache Even Sequence. Are there terms in this sequence which are 2*p 
where p is a prime (or pseudo prime). With a small change in the program 
used for the Smarandache Odd Sequence we can easily find out Strangely 
enough not a angle term was found to be of the form 2p. 

V.4 The Smarandache Prime Sequence 

The Smarandache Prime Sequence is generated by {2,3,5,7,11, Again 
we ask: - How many are primes? - and again we apply the method of finding 
the number of primes/pseudo primes among the first 200 terms. 

There are only 4 cases to consider Terms #2 and #4 are primes, namely 23 
and 2357. The other two cases are: term #128 which is a 355 digit number 
and term #174 which is a 499 digit number. 

#128 

23571 113171 9232931 3741 4347535961 6771 73798389971 01 1 031 071 09 1 1 31 271 31 1 371 39 
1 49 1 5 1 1 571 631 671 731 791 81 1 9 11 931 971 9921 1 223227229233239241 25 1 257263269271 2 
77281 283293X731 131 XI 7X1 X734734935335936737X79383389397401 40941 9421 43 



‘Determine length of U% 
‘Add on U% 

‘Determine E% (=n) 
‘Determine A (=a) 

Check if a is a n/h power 
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1 43343944344945746 1 463467479487491 499503509521 523541547557563569571 577587 
593599601607613617619631641643647653659661673677683691701709719 

Table 1 c. The first 20 terms of the Smarandache Prime Sequence 



2 

23 
235 
2357 
235711 
23571113 
2357111317 
235711131719 
23571113171923 
2357111317192329 
235711131719232931 
23571113171923293137 
2357111317192329313741 
235711131719232931374143 
23571 1 1 3 1 71923293 1 374 1 4347 
23571 1 1317192329313741434753 
23571 1 131 719232931374143475359 
23571113171923293137414347535961 
2357111317192329313741434753596167 
2357 1 1 1 3 1 71923293 13741 434753596 1 6771 
235711131719232931374143475359616771 73 



#174 

2357111317192329313741434753596167717379838997101103107109113127131137139 
1491511571631671731791811911931971992112232272292332392412512572632692712 
77281 28329330731 131331 7331337347349353359367373379383389397401 40941 9421 43 
1 433439443449457461 463467479487491 499503509521 523541547557563569571 577587 
593599601 60761361 7619631 641 643647653659661 673677683691 701 70971 97277337397 
43751 757761 76977378779780981 182182382782983985385785986387788188388790791 
1919929937941947953967971977983991997100910131019102110311033 



Are these two numbers prime numbers? 
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Chapter VI 
Oq the Harmonic Series 



VL1 Comparison of a few seqoacei 

t>v harmonic series 1 ~+ — is so important and so intensely 

studied that nothing mens can be said about it - or can there? We may always 
pose a few questions. 

® 1 

We know that the series is divergent Furthermore the sates Z — is 

convergent for a>l and divergent for a£l. Since it is the “borderline” it may 
be interesting to examine the following questions: 



Question 1. what is the smallest number of terms m we need in order to 
tnairp the qim greater than a given positive integer w? In other words the 

m 1 

smallest m for which £ r- w - 

k=i* 



Question 2. What is the smallest m for which £— r^w? 

k=i k* 

Using the notations for the Smarandache function and the square residues 
intirwhwfd in chapter H.2,3 we pose the correspo n di ng questions: 

m 1 

Question 3. What is the smallest m for which Z w? 

k=tS (k) 

m l 

Question 4. What is the smallest m for which I , £w? 

k=iS r ‘(k) 



Table 1 shows the values of m required in each case for w=l, 2, 3, ... 12. 
diagram i ln(m) is plotted against w. A simplified version of this diagram 

® 1 

shown cm the cover. We see that the series I — is very “sensitive” to 

k=lk 
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changes in the value of a. 



Z — - — 

k=.S r 2 (k) 



is on the convergent side of the graph for 



m 1 

1/k and Z — rr : on the divergent side. 

k=iS 2 (k) 



Table 1 . Comparison of sequences. 



w 


Series 1: 
1/k 

Values of ro 


Series 2: 
l/S 2 (k) 

Values of m 


Series 3: 
1/S r 2 (k) 

Values of m 


Series 4: 
l/k 1 - 02 

Values of m 


1 


2 


24 


2 


2 


2 


4 


144 


9 


4 


3 


11 


462 


54 


12 


4 


31 


1045 


243 


35 


5 


83 


1944 


729 


102 


6 


227 


3200 


2048 


311 


7 


616 


4862 


6561 


966 


8 


1674 


6912 


16384 


3082 


9 


4550 


9477 


37746 


10115 


10 


12367 


12480 


93312 


34167 


11 


33617 


16008 


209952 


118971 


12 


91380 


20020 


497664 


427698 




Diagram 1. The behaviour of some sequences: ln(m) plotted against w. 
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VL2 Integers represented as sums of terras of the harmonic series 

Given a positive integer w. Can we represent w as a sum of a finite number of 
different terms chosen from the harmonic series? For w=l we have the simple 

answer l , but what about w=2, w=3 etc. 

2 3 6 

Method 1: k-perfect numbers. By definition a k -perfect number is a number 
n which satisfies the equation cr(n)=k-n, where o(n) is the sum of divisors of 

(1 and n included), k = therefore gives the desired representation of the 
n 

integer k. The result of a computer search for k-perfect numbers is given in 
table 2. A column is included for the number of divisors m since this is equal 
to the number of terms required for the representation. The numbers are given 
in factor form as this is more interesting than the decimal representation. A 2- 
perfect number is in the literature simply referred to a perfect number. They 
have attracted a lot of interest - The problem of perfect numbers* a favorite 
with ancient Greeks ; owes its origin to the number mysticism of the 
Pythagoreans (quote from Elementary Number Theory, Uspensky Headet). 



Table 2. k-perfect numbers a m = the number of dtasors. 



k 


m 


n 


2 


4 


2*3=6 


2 


6 


2 2 -7=28 


2 


10 


2*-31=496 


2 


14 


2*127=8128 


2 


26 


2 12 -8191=33550336 


2 


34 


2 1S -131071=8589869056 


2 


38 


2 18 -5242 87=1374 38 69132 8 


2 


62 


2 30 -2 14 74 83 64 7=2305 843008 139952128 


~3 


16 


2 3 -3-5=120 


3 


24 


2 s -3-7=672 


3 


288 


2 8 -5-7-l 9-37-73=45 98 182 40 


3 


80 


2 9 -3-l 1-31=523776 


3 


224 


2 13 -3-l 1-43-127=1 476304 896 


3 


480 


2 14 -5-7-l 9-31-15 1=5 1001180160 
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Table 2. continued. 



k 


M 


n 


4 


216 


2 2 -3 2 -5-7 2 - 13-1 9=2 1785 40 


4 


96 


2 3 -3 2 -5-7-13=32760 


4 


96 


2 s 3 3 -5-7=30240 


4 


384 


2 7 -3 3 -5 2 -17-31=45532800 


4 


320 


2 9 -3 3 -5-ll-31=23569920 


4 


480 


2 9 -3 2 -7- 11-1 3-3 1=1 4 2 9908 4 8 


4 


1056 


2 10 -3 3 -5 2 -23-31-89=43861478400 


4 


9984 


2 25 -3 3 -5 2 -l 9-3 1-683-273 1-8 191 


5 


1920 


2 7 -3 4 -5-7-ll 2 - 17- 19=14182439040 


5 


2304 


2 7 -3 5 -5-7 2 -13-17-l 9=3 19983955 20 


5 


5280 


2 10 -3 4 -5-7-ll 2 -l 9-23-8 9 


5 


6336 


2 10 -3 5 -5-7 2 -13-19-23-89 


5 


3456 


2 11 -3 3 -5 2 -7 2 -13-19-31 


5 


9216 


2 11 -3 s -5 3 -7 3 -13 3 -17 


5 


43008 


2 20 -3 3 -5 3 -7 3 -13 3 -17- 127-3 3 7 


5 


709632 


2 21 -3 6 -5 2 -7-19-23 2 -31-79-89-137-547-683-1093 


5 


94208 


2 22 -3 7 -5-7-ll-l 9-41-4 7- 151-1 97-178481 


6 


1658880 


2 17 -3 4 -5 3 -7 3 -ll 2 -13 2 -19 3 -31-37-61-73-i81 



For w=3 we have a representation by 16 terms: 

,.1111111 1 1 1 1 1 1 l l 

3= 1+— + — + — H 1 1 1 1 1 1 H f- 1 1 -( 

2 3 4 5 6 8 10 12 15 20 24 30 40 60 120 

The smallest number of terms required for w=5 is 1920 and for 6 a staggering 
1,658,880 terms. For w=7 we have no representation at all. We would first 
have to find 7-perfect number. 



No great results. It’s like having used a sledge hammer to kill a mosquito and 
missed the target. But k-perfect number are interesting in themselves. 



Method 2: Trivial Expansions. A term in the harmonic series can always 

be replaced by two other different terms from the harmonic series: 
1 k+1 1 ... 

k = k(k + l) = ~k + l + k(k + \) ‘ GlVCn 3X1 ex P ansion of m mte S er w into terms 
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of the harmonic series it is possible to construct infinitely many such trivial 
expansions of w. If however, in a given expansion of w, each term is rep l aced 
by two other terms as above and furthermore the re pl ace ment process is 
carried out until all terms in the new represe n tation are different from one 
another and also from the terms in the original r e pr ese ntat ion then a new 
representation has been obtained. Obviously such repres ent a tion s for Wi and 
W 2 with no common terms can be added together to form a representation of 
W=Wi+W2- 



Example: 

i 1 1 1 

1= 2 + 3 + 6 










Expand: 


k 


-» 


k+l 


k(k+l) 


2 




3 


6 


Not to be used (repeated) 


3 


—> 


4 


12 




6 




7 


42 





3 and 6 are repeated and since their expansion gives 4, 12 and 7, 42 we have 
to continue the expansion: 



4 


— ► 


5 


20 


12 


— ► 


13 


156 


7 


-* 


8 


56 


42 




43 


1806 



From this we can write the following representation of 2: 



2=i+I+I+I 

4 5 7 



11111111 1 
— + -+■ — + ■f " — + — — Hh *~~ 4 — + 

8 12 13 20 42 43 56 156 1806 



Adding our original representation of 1 we get the following representation 3: 
1 1 1 1 1 1 1 111111 11 
3 = 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 12 13 20 42 43 56 156 1806 



As in the case of 3 -perfect numbers we have a representation of 3 with 16 
terms. The minimum number of terms for which a repres entatio n is possible 
is 11 as can be seen from table 1. Let’s try to get closer to this. 



Method 3, M aximum density. For a given value of w table 1 gives the 
minimum number of terms from the harmonic series required to form a 
representation For a given value of w=3 this number is m. We can therefore 
try to choose the first m-1 terms as dose together as possible, i.e. we farm 
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s=w- Z- 
k=l* 



where s is a small fraction which we have to expand in harmonic terms. As 
before we apply the idea to w=3. We get 



^llllllllll 1 1 

3 = 1 +—+—+—+-+— + —+—+—+ — + — + + 

23456789 10 15 230 57960 

Bravo! Only 13 terms. We had 16 with the other two methods. But no 
pleasure lasts forever. The nasty little fraction s will give no end to problem 
when we try large values for w. Although small s will have large numerators 
and denominators as we shall soon see. Maybe better get the sledge hammer 
out again and go in search for the 7-perfect number. 



VL3 Partial sums of the harmonic series as rational numbers. 



n 1 

We consider consecutive partial sums S n = £ — expressed as rational 

k=l k 

numbers reduced to lowest terms. 



Question 1. Can there be a repetition of numerators? 

Assume S n _! and S n =-^- where (a,b)=l and (a,bi)=l. We then have 

b bj 

— which is equivalent to an(b-b 1 )» but (a,bbi)=l and the latter 

bj b n 

equation is therefore impossible and the assumption is wrong. We conclude: 



Two consecutive sums Sn_i and Sn cannot have equal numerators when 
reduced to lowest terms. 



A & 

A necessary condition for fidig two sums S n = — and S m = — — , m>n+l with 

b n b m 

equal numerators a n =a m when in lowest terms is therefore that ri\ exists so that 
n<ni£m and a nj <a ni _! . Such numerators exist Examples: 

a € =49 

a 7 =137 

ai8=14274301 

an=4214223 
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3*9=360968703235711654233892612988250163157207 
aioo=l 4 4 66 63 62795 2035 11 602215 18 043 104 1314 47711 
etc 

Now, let p be a prime then because Sn is a monotonously increasing 

function of n and 

s *P-1 , l a P-lP + b P-l 

p bp_, p bp_,p 

cannot be reduced to lower toms. 



Conjecture: Let pj and p^-i denote consecutive primes. The numerators in die 
interval p^n^p^.,, after reducing to lowest terns, are all greater than » p ^ -1 , 

i.e. when a term of the type 1 is added the numerator takes a leap to a higher 

P 

value from which, although it may decrease, it will never drop to a value 
below the value before the leap. 



Question 2. What is the most frequently occurring denominator in the 

sequence S. = 1+— +— + ... +— when reduced to lowest terms? 

0 2 3 n 



Let S a _, =— where (a,b)=l. Thai S„=^-+— ; S„ = ^ ~ ^. Consider the 
b b n on 

greatest common divisor of b and n, d=(b,n). Then it is seen immediately that 

the fraction can be reduced by d. Pitt b=bid and n=nid and consequattly 

ani+b, 

^ b,n,d 

If a further reduction is possible then it must be by a factor in d or by d itself 
because (nj,bi)=l and (bi,a)=l. 



If n=p is a prime number then (b,n)=l and the denominator os Sn will have p 
as a factor. The denominator of S* is therefore unequal to all previous 
denominators and is at most the first in a series of equal denominators. Sums 
with equal denominators must therefore be contained between 

1 . _ .11 1 



o , 1 1 

s », 3+- 



“■ v.= 1 + rr~ V 



86 




where ft and p^i are consecutive primes. In table 3 the number of equal 
denominators is denoted q and the corresponding number of digits of the 
denominator L. 



Table 3. Equal denominators for n<1000. 



Pi 


Pi+i-1 


q 


L 


89 


96 


8 


39 


139 


148 


10 


62 


317 


330 


14 


140 


891 


906 


16 


389 



I hope this bo ok has been more meaningful than reading the number below 
which is the 389-digit denominator which occurs 16 times in consecutive 
partial sums of the harmonic series for n<1000. 



4179982336319706196962068134998944512773598562712 

1827499073959690180767635682039557104883261454628 

9717289725066739371544332075631263681868936793360 

4957159732797690125645682733080622438318699474351 

6944805893417522011858358480060385932366764886574 

4940462245050719348300672374671465097882430105554 

5638386382683543980298959618721940328647837326571 

6039813081070459435310984219943740583597120000 





